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ABSTRACT 

The  present  work  is  concerned  with  a  theoretical  and  an  experi¬ 
mental  treatment  of  non- equilibrium  flow  of  a  chemically  reacting  gas 
mixture  through  a  shock  wave.  The  two  chemical  reactions  given  by 
N2O4  +  M  ■  2NC>2  +  M  =  2NO  +  O2  +  M,  where  M  represents  the  inert 
carrier  gases  argon  or  nitrogen,  were  studied.  The  experiments  were 
carried  out  in  a  temperature  controlled  shock  tube  and  a  light  absorp¬ 
tion  technique  permitted  the  time  dependent  concentration  of  the  species 
NO2  to  be  determined.  For  shock  strengths  where  the  temperature 
did  not  exceed  400°K  only  the  first  chemical  reaction  took  place. 
Stronger  shock  waves  excited  both  chemical  reactions  with  tempera¬ 
tures  up  to  2100°K,  and  for  this  reason,  flow  fields  with  two  non¬ 
equilibrium  modes  could  be  investigated.  Since  the  relaxation  times 
of  these  two  reactions  were  different  by  about  three  orders  of  magni¬ 
tude,  they  could  be  experimentally  uncoupled.  A  study  of  the  reaction 
mechanisms  and  rate  constants  for  both  reactions  was  carried  out. 

In  general,  the  experimental  work  agreed  well  with  an  analytical 
description  of  the  non- equilibrium  flow  fields,  and  with  the  work  of 
previous  investigators  at  conditions  where  direct  comparisons  could 
be  made.  At  shock  strengths  exciting  only  the  first  chemical  reaction 
the  complete  picture  of  a  non-equilibrium  flow  field  with  only  one  non¬ 
equilibrium  mode  could  be  investigated.  In  this  situation,  the  shock 
strengths  were  varied  from  weak,  fully  dispersed  waves  to  strong, 
partly  dispersed  waves.  It  is  believed  that  fully  dispersed  waves  were 
observed  for  the  first  time  in  a  chemically  reacting  system,  in  con¬ 
trast  to  previous  observations  involving  vibrational  non- equilibrium. 
Finally,  it  is  demonstrated  that  fully  dispersed  waves,  with  their 
inherently  long  relaxation  lengths  and  small  gradients  in  state  vari¬ 
ables,  may  be  used  as  a  particularly  simple  technique  for  the  chemical 
kinetic  study  of  flows  with  one  non- equilibrium  mode. 
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SECTION  I 
INTRODUCTION 


At  the  temperatures  achieved  in  hypersonic  flew  fields,  processes 
such  as  vibrational  excitation,  thermal  dissociation,  chemical  reactions, 
ionization,  radiation,  and  electronic  excitation  may  occur.  In  air  for 
example,  at  Mach  numbers  above  5,  Stupochenko  (1)  lists  six  different 
chemical  reactions  alone.  For  the  short  flew  times  ordinarily  character¬ 
istic  of  high  velocity  gas  dynamics,  and  in  particular  at  higher  altitudes, 
the  assumption  that  the  fluid  is  everywhere  in  local  thermodynamic  equili¬ 
brium  is  not  always  valid,  and  the  non-equilibrium  mechanisms  must  be  in¬ 
cluded  in  a  description  of  the  flow  field. 

The  large  number  of  simultaneous  reactions  in  air,  and  the  high 
temperatures  required  to  excite  these  reactions  causes  experimental  in¬ 
vestigations  in  air  to  be  complicated.  In  order  to  investigate  certain 
aspects  of  the  fundamentals  of  non-equilibrium  flow  in  the  laboratory,  a 
different  gas  system  may  be  chosen  in  which  the  temperatures  necessary  to 
induce  excitation  are  readily  obtainable,  and  in  which  the  non-equilibrium 
mechanisms  under  study  may  be  resolved  and  uncoupled  from  each  other.  For 
these  reasons  (2)  a  flew  was  chosen  in  which  the  following  reactions  occurs 


N204  +  M^  2N02  +  M  *  2N0  +  02  +  M 


(1) 


where  M  represents  an  inert  gas,  either  nitrogen  or  argon.  Also,  the 
reactant  mole  fraction  did  not  exceed  10%.  The  flow  fields  were  produced 
in  a  temperature  controlled  shock  tube  in  which  shock  waves  of  different 
strengths  were  propagated  through  a  knewn  equilibrium  mixture  of  the  gases 
of  Equation  (1) . 

The  work  was  concerned  with  three  major  objectives. 

(A)  One  aim  was  to  investigate  the  complete  picture  of  a  non-equilibrium 
flew  field  for  shock  waves  varying  from  weak,  fully  dispersed  waves 
to  strong,  partly  dispersed  ones.  The  waves  were  to  be  propagated 
through  a  gas  mixture  with  a  single  active  non-equilibrium  mode  whose 
reaction  mechanism  and  rate  constants  are  well  understood. 

(B)  A  second  objective  was  to  include  a  flow  field  with  two  non-equilibrium 
processes.  Since  the  relaxation  rate  for  the  second  reaction  of  Equa¬ 
tion  (1)  is  much  slower  than  that  for  the  first  reaction,  these  two 
processes  could  be  experimentally  separated. 

(C)  The  third  goal  was  to  repeat  and  extend  measurements  of  the  rate  con¬ 
stants  for  both  reactions.  In  particular  an  extension  of  results  at 
different  temperatures  and  reactant  concentrations  than  currently 
available  in  the  literature  appeared  to  be  important.  Existing  rate 
equations  were  to  be  examined  in  this  process  and  compared  with  the 
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results  of  References  3  through  9. 

A  shock  wave  may  be  defined  as  a  non-equilibrium  layer  separating 
two  equilibrium  states  in  a  gas.  As  the  gas  flows  through  this  non¬ 
equilibrium  layer,  molecular  collisions  change  the  state  of  the  gas  from 
one  equilibrium  state  to  the  second.  Each  degree  of  freedom  of  a  mole¬ 
cule  requires  a  certain  number  of  collisions  to  achieve  the  new  equili¬ 
brium  state.  The  time  required  for  these  collisions  to  take  place  is 
kncwn  as  the  relaxation  time  and  each  degree  of  freedom  has  its  own 
characteristic  relaxation  time. 

A  characteristic  time  may  also  be  ascribed  to  the  flow.  For 
example,  in  our  flews  the  time  required  for  a  particle  to  traverse  the 
non-equilibrium  layer  of  a  shock  wave  may  be  chosen.  The  ratio  of  the 
relaxation  time  to  this  particle  flew  time  has  often  been  called  the 
Damkohler  parameter  D  .*  For  D  e  <*>  ,  the  process  is  frozen  since  there 
is  no  change  in  the  degree  of  freedom  across  the  region  of  interest.  For 
D  =  0  ,  the  process  is  in  equilibrium  since  the  degree  of  freedom  relaxes 
instantaneously.  For  most  applications  in  flew  systems,  conditions  at 
these  two  limits  provide  an  adequate  description  of  the  flow  fields.  How¬ 
ever,  for  situations  where  D  =  1  ,  non-equilibrium  effects  must  be  con¬ 
sidered.  For  the  system  under  consideration  in  this  work, 


t 


T' 


V 


<< 


T 

a 


<<  Te 


(2) 


in  which  the  subscripts  refer  to  the  translational,  rotational,  vibrational 
and  the  two  chemical  relaxation  processes  respectively.  The  translational 
and  rotational  modes  are  here  listed  together  since  they  are  close  to  each 
other  with  respect  to  the  relaxation-  time  for  vibrational  adjustment. 

In  the  region  of  interest  for  the  first  chemical  reaction,  the 
translational,  rotational  and  vibrational  degrees  of  freedom  may  be  con¬ 
sidered  in  equilibrium.  In  addition,  the  second  reaction  is  frozen.  In 
the  region  of  interest  for  the  second  reaction  at  higher  temperatures,  all 
the  other  degrees  of  freedom,  including  the  first  reaction  may  be  considered 
in  equilibrium. 

An  essential  feature  of  a  reacting  gas  mixture  is  the  fact  that 
the  mixture  behaves  as  a  dispersive  medium  from  the  viewpoint  of  acoustics. 
There  is  no  unique  acoustic  velocity  and  frequency  enters  as  a  parameter 
(11) .  The  acoustic  equations  for  flew  with  vibrational  or  chemical  non¬ 
equilibrium  are  developed  by  Vincent!  and  Kruger  (12) .  Two  limiting 
sound  speeds  may  be  defined  corresponding  to  the  frozen  and  equilibrium 
limits  (13) .  The  frozen  sound  speed  af  propagates  as  if  all  processes 
in  the  gas  were  frozen.  The  equilibrium  sound  speed  a@  which  is  always 
less  than  af  ,  (14) ,  propagates  as  if  all  processes  in  the  gas  were  in 


*Wegener  and  Buzyna  (10)  give  seme  of  the  history  of  this  parameter  and 
its  significance  in  the  description  of  non-equilibrium  flow  fields. 
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equilibrium. 

These  two  limiting  sound  speeds  (12)  are  expressed  as  follows . 
The  frozen  sound  speed  a^  is  given  by 


4ft 


s,a,8 


(3) 


where  a^2  =  yKT  with  y  =  Cp/C^  if  the  species  are  individually  ther¬ 
mally  perfect.  Furthermore,  thevequilibrium  sound  speed  aQ  can  be 
shown  to  be 


4ft 


-  +  h 
P 


l  a 
a  p 


(4) 


Here,  the  asterisk  refers  to  the  equilibrium  values  of  the  quantities, 
and  the  variables  with  subscripts  indicate  partial  derivatives.  These  de¬ 
rivatives  were  first  evaluated  by  Wegener  and  Cole  (15) ,  specifically  for 
the  first  reaction  of  Equation  (1) .  As  mentioned  first  by  Be  the  and 
Teller  (16)  the  possibility  of  shock  waves  traveling  between  these  two 
limiting  sound  speeds  exists/  and  Griffith  and  Kenny  (17)  demonstrated  the 
existence  of  these  waves  for  vibrational  relaxation  in  002. 

In  this  work,  these  waves  were  observed  for  the  first  time  in  a 
chemically  reacting  system.  Waves  traveling  between  the  limiting  sound 
speeds  were  called  fully  dispersed  by  Lighthill  (18) ,  while  waves  propa¬ 
gating  at  velocities  above  the  frozen  sound  speed  were  called  partly  dis¬ 
persed. 


In  a  partly  dispersed  wave,  there  is  a  normal  shock  jump  in  which 
the  rapidly  adjusted  degrees  of  freedom  relax.  In  our  case  they  are  trans¬ 
lation,  rotation,  and  vibration.  In  this  short  interval  the  slow  degrees, 
i.e.  the  chemical  reaction,  remains  frozen.  After  the  normal  shock  jump, 
the  slew  degrees  of  freedom  proceed  to  the  final  equilibrium  conditions. 

As  the  shock  strength  is  decreased,  the  overall  change  in  state  caused  by 
both  the  normal  shock  junp  and  the  chemical  reaction  decreases.  However, 
the  chemical  reaction  produces  an  increasing  fraction  of  the  total  change 
in  state.  Finally,  for  a  shock  wave  propagating  exactly  at  the  frozen 
sound  speed,  the  normal  shock  jump  goes  to  zero,  and  the  entire  change  in 
state  is  due  solely  to  chemical  relaxation.  In  other  words,  if  the  flew 
velocity  is  related  to  the  frozen  sound  speed,  at  a  Mach  number  of  one 
(M  =  U/af)  the  normal  shock  vanishes.  Belcw  the  frozen  sound  speed,  fully 
dispersed  waves  provide  a  continuous  transition  between  the  two  equilibrium 
states. 
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Fig.  1  Schematic  Pressure  Change  for  Partly  and  Fully  Dispersed  Shock  Waves; 
Flow  Direction  from  Left  to  R  ight 


Figure  1  shews  the  typical  behavior  for  pressure  of  partly  and 
fully  dispersed  shock  waves  for  an  endothermic*  reaction.  The  partly  dis¬ 
persed  wave  is  characterized  by  the  normal  shock  jump  followed  by  the  re¬ 
gion  of  non-equilibrium  relaxation.  The  fully  dispersed  wave  exhibits  a 
continuous  non-equilibrium  transition  between  the  two  equilibrium  states, 
with  no  normal  shock  junp  apparent. 

There  are  several  reasons  why  the  use  of  fully  dispersed  waves  is 
ideally  suited  for  a  study  of  chemical  kinetics.  The  chief  reason  is  that 
the  temperature  of  the  reaction  can  be  very  accurately  determined  since  the 
initial  temperature  of  the  wave  is  measured  statically.  Typically  in  our 
work,  the  total  temperature  change  across  the  wave  is  only  about  ten  de¬ 
grees  Kelvin.  This  accuracy  in  temperature  is  important  since  rate  con¬ 
stants  are  temperature  sensitive  since  they  obey  laws  such  as  k  =  A  exp 
(-  B/T)  where  A  and  B  are  constants.  Another  reason  is  that  the 
smooth  transition  between  states  in  these  waves  facilitates  the  measure¬ 
ments  of  slopes.  Also,  the  lew  velocities  and  long  relaxation  times  in¬ 
crease  the  resolution  of  the  experiment.  Finally,  the  small  gradients  in 
the  flew  variables  make  these  waves  amenable  to  simplified  theories. 

The  dispersion  owing  to  the  chemical  reaction  proceeding  "slowly" 
tends  to  spread  the  shock  wave  out  in  time  (or  space) .  This  effect  is 
balanced  by  the  convective  terms  in  the  conservation  equations  which  tend 
to  steepen  the  wave.  In  the  steady  state,  the  particles  flowing  into  the 
shock  front  pass  through  the  same  sequence  of  states  as  the  preceding  ones. 
The  distribution  of  the  variables  across  the  shock  front  thus  forms  a  fixed 
picture  Which  moves  during  this  period  as  an  entity  together  with  the  front 
(14) . 


*Heat  is  absorbed  by  the  chemical  reaction  as  typically  in  dissociation. 
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SECTION  II 

EQUILIBRIUM  PROPERTIES  OF  THE  NITROGEN  DIOXIDE 
SYSTEM  AND  EQUATIONS  OF  STATE 


The  chemical  reactions  observed  in  the  non-equilibrium  flew  field 

were: 


N204  +  M  t  2N02  +  M  t  2N0  +  02  +  M  (5) 


in  which  M  ,  an  inert  component,  was  either  nitrogen  or  argon. 
The  mole  fractions  of  the  constituents  are;* 


XoU  -  a) 

^204  =  1  +  aX  +  BX 
o  o 

2XQ(a  -  B) 

^2  =  1  +  aX  +  BX 
*  o  o 

2Xq& 

^0  =  1  +  aX  +  BX 
o  o 


*0,  1  +  aX  +  BX 

o  o 

and, 

1  -  X0 

=  1  +  aX  +  BX 
o  O 


In  terms  of  these  quantities,  the  laws  of  mass  action  for  the  reactions 
are: 


*In  all  the  equations  to  be  given,  values  of  numerical  constants  etc.  are 
consistent  with  the  units  given  in  the  list  of  symbols. 
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K  -p  X°  4(g  -  8) 2 

pi  F  1  +  ctXQ  +  0XQ  1  -  a 

and 

X 

K  =  n  O  3 

p2  P  1  +  <xXo  +  0Xo  (a  -  3)2 


(7) 


(8) 


vdiere  subscripts  1  arid  2  refer  to  the  first  and  second  parts  of  the  reac¬ 
tion  respectively. 

Kpl  has  been  measured  by  Bodenstein  &  Boes  (6) .  The  results  were 
compiled  by  Giauque  and  Kemp  (19)  and  by  Strehlow  (20) .  The  latter’s 
critical  evaluation  yields. 


iogio  Kpl  =  -  +  9.187  (9) 


where  T  is  in  degrees  Kelvin  and  the  equilibrium  constant,  Kpi  in 
atmospheres  pressure.*  This  part  of  the  reaction  takes  place  in  the  range 
200  °K  <  T  <  400  °K. 

The  second  part  of  the  reaction  of  Equation  (5)  was  investigated 
by  Bodenstein  and  Lindner  (6)  and  also  compiled  by  Giauque  &  Kenp  (19) . 
Fran  their  data  we  find  in  the  units  of  Equation  (9) : 


logio  Kp2  =  -  ^  +  7.803  (10) 


This  reaction  takes  place  in  the  range  500  °K  <  T  <  1200  °K.  Much  further 
information  including  thermodynamic  and  transport  properties  is  reviewed  in 
References  21,  22  and  23. 

At  the  tenperature  and  relatively  low  pressures  used  in  the  shock 
tube  in  the  following  experiments,  the  constituents  of  the  gas  may  indivi¬ 
dually  be  considered  as  thermally  perfect  gases  since  intermolecular  forces 
are  negligible. 

If  W  ,  the  molecular  weight  of  the  mixture  is  considered  as  a 
variable  of  the  non-equilibrium  flow  field,  the  thermal  equation  of  state 
for  the  mixture  becomes 


p  _  RT 
p  “  W 


(11) 


*1  atmosphere  =  760  torr. 
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where  R  is  the  universal  gas  constant.  We  have:* 


M  =  WA  +  92’02  V  +  46-01  xro2  +  30-01  +  32'00  ^  (L2> 


The  caloric  equation  of  state  is 


otX. 


h  •  W  'CP(T  -  V  +  13846  X  +  aX  +  BX  + 


27020 


p '  o'  1  +  ctX  +  $X  “ ' 1  +  oX  +  $X 

c  o  o  o  < 


-]  (13) 


where 


Cp  CpM  +  ^NjOi*  CpN204  +  XN02  CpN02  +  ^  CpN0  +  X02  Cp02  (14) 


Values  of  the  specific  heats  of  the  constituents  and  of  the  heats 
of  reaction  were  taken  from  References  24  and  22  respectively. 


*With  the  tabulated  values  of  molecular  weights  of  the  National  Bureau 
of  Standards. 
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SECTION  III 

ANALYSIS  OF  ONE-DIMENSIONAL  NON-EQUILIBRIUM  FLOW 
THROUGH  A  SHOCK  WAVE 


The  shock  tube  (25)  is  a  length  of  tubing  divided  into  two  sec¬ 
tions  by  a  breakable  diaphragm  as  shown  in  Figure  1.1  in  Appendix  I,  where 
initially  a  pressure  difference  exists  between  the  sections. 

Figure  2  is  a  characteristic  x-t  diagram  of  the  flow  in  a  shock 
tube.  When  the  diaphragm  is  broken,  the  high  pressure  driver  gas  expands 
into  the  driven  tube.  The  front  of  the  driver  gas,  called  the  contact  sur¬ 
face,  acts  as  a  piston  and  drives  a  shock  wave  in  front  of  it.  A  rare¬ 
faction  wave  meanwhile  passes  through  the  driver  gas  and  lowers  its  pres¬ 
sure  and  temperature.  The  testing  time  available  at  the  observation  window 
depends  on  the  lengths  of  the  two  tube  sections.  It  is  limited  by  either 
the  passage  of  the  contact  surface,  the  passage  of  the  reflected  shock  wave, 
or  the  passage  of  the  reflected  rarefaction  wave. 

In  a  shock  tube,  at  a  sufficient  distance  downstream  of  the  dia¬ 
phragm,  the  shock  wave,  across  which  the  variables  of  state  jump  and  re¬ 
main  fixed  relative  to  each  other,  propagates  at  a  constant  velocity  into 
the  quiescent  gas  in  front  of  it.* 


Fig.  2  Characteristic  Diagram  of  Non-Equilibrium  Flow  in  a  Shock  Tube 


In  Figure  2,  state  1  applies  to  the  initial  equilibrium  gas  mix¬ 
ture  at  rest.  State  2  occurs  after  the  normal  shock  jump  in  which  the 
translational ,  rotational,  and  vibrational  degrees  of  freedom  relax  and 
the  chemical  reactions  are  frozen.  At  state  3  the  first  reaction  of  Equa¬ 
tion  (1)  reaches  equilibrium  while  the  second  reaction  is  still  frozen. 
Finally,  at  state  4  all  degrees  of  freedom  are  again  in  equilibrium.  Typi¬ 
cal  values  of  our  relaxation  times  are  tt  ,  xR  *  10“  3  ps  ,  tv  =  0.3  js  , 


*For  the  actual  starting  process  in  a  shock  tube  see  Spence  (26) . 
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and  =  10  us  at  T  =  320  ®K  and  p  =  250  torr. 

In  the  experiments  considered  in  this  work,  a  temperature  range 
of  280  <  T  <  2100  °K  and  a  pressure  range  of  200  <  p  <  1000  torr  was  in¬ 
vestigated.  Since  the  gradients  in  the  flew  properties  were  not  great, 
except  in  the  normal  shock  jump,  the  effects  of  viscosity,  heat  conduction, 
and  diffusion  wejre  neglected.  Also,  a  one  dimensional  flow  approach  was 
valid  since  boundary  layer  effects  that  slow  dewn  the  shock  were  shown  to 
be  negligible. 

In  the  analysis  that  follows  (36)  a  coordinate  transformation  is 
made  in  which  the  shock  is  fixed  while  the  flow  entering  the  shock  front 
has  the  measured  shock  velocity  Uj  . 

The  steady  state,  one  dimensional  conservation  equations  are. 


Mass: 

PjUi 

p 

Q. 

II 

(15) 

Momentum: 

Pi  + 

~cu 

ii 

CM 

D 

H 

Q. 

+  p'U'2 

(16) 

Energy: 

hi  + 

|U!2  =  h' 

+  IU'2' 

(17) 

The  primes  denote  conditions  anywhere  in  the  flow  field.  For  given  initial 
and  boundary  conditions,  these  three  equations  are  not  sufficient  to  cal¬ 
culate  the  four  variables  p  ,  U  ,  p  ,  and  h  ,  and  an  equation  of  state 
is  necessary  to  relate  these  variables.  If  a  gas  is  in  unconstrained 
thermodynamic  equilibrium,  two  variables  are  sufficient  to  describe  the 
thermodynamic  state  of  the  medium.  For  example, 


h*  =  h*  (p,  p)  (18) 


If  the  gas  is  not  in  equilibrium,  then  additional  variables  of 
state  are  required  for  each  non-equilibrium  mode.  The  canonical  equation 
of  state  (12)  beccmes, 


h  =  h  (p,  p,  0^,  a,  B)  (19) 


Here  Cp,  the  specific  heat,  refers  to  the  translational,  rotational  and  vi¬ 
brational  modes,  and  a  and  8  refer  to  the  two  chemical  nan-equilibrium  de¬ 
grees  of  freedom.  At  equilibrium.  Equation  (19)  reduces  to  Equation  (18) ' 
since  the  non-equilibrium  variables  became  functions  of  p  and  p  alone.  The 
caloric  and  thermal  equations  of  state  for  the  gas  system  under  consideration 
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are  qiven  by  Equations  (13)  and  (11) . 

These  equations  together  with  the  conservation  equations  and  for 
the  given  initial  conditions  of  Pi  t  Pi  /  Ui  ,  and  hi  are  sufficient  to 
solve  the  flow  conditions  at  the  states  2,  3,  and  4  of  Figure  2. 

State  2  occurs  after  the  normal  shock  junp  in  which  the  transla¬ 
tional,  rotational,  and  vibrational  degrees  of  freedom  adjust.  As  a  first 
approximation  to  this  state,  Cp  ,  a  ,  and  8  are  frozen  and  the  caloric 
equation  of  state  for  a  thermally  perfect  gas  becomes  simply. 


h' 


=  —  C  T' 
W  pi  1 


(20) 


where  the  reference  enthalpy,  hx  ,  has  been  set  equal  to  zero.  The  solu¬ 
tion  of  Equations  (15) ,  (16) ,  (17) ,  (11)  and  (20)  represents  the  well- 
known  Rankine-Hugoniot  solution  (27)  for  a  thermally  and  calorically  per¬ 
fect  gas,  and  one  obtains  p'  ,  p'  ,  U'  ,  and  h'  .  We  define. 


C  dT  (21) 

P 

Tl 

where  Cp  is  the  relaxed  heat  capacity  at  T  .  With  this,  the  caloric 
equation  of  state  becomes, 


h2  =  i  Cp  (T2  -  TX)  (22) 


Setting  h2  equal  to  h1  initially.  Equations  (15) ,  (16) ,  (17) ,  (11)  and 
(22)  can  be  solved  by  iteration  to  obtain  the  conditions  for  state  2. 

At  state  3,  the  first  reaction  of  Equation  (1)  goes  to  equilibrium 
while  the  second  reaction  is  considered  to  remain  frozen.  To  obtain  con¬ 
ditions  at  this  point  the  law  of  mass  action.  Equation  (7) ,  is  employed. 
From  Equations  (7)  and  (9) ,  the  degree  of  dissociation  is  calculated  for 
particular  values  of  pressure  and  temperature .  Initial  values  of  state  3 
are  assumed,  and  with  Equation  (7)  the  Equations  (15) ,  (16) ,  (17) ,  (11)  and 
(13)  may  again  be  solved  by  iteration  to  obtain  conditions  for  state  .3. 

The  -final  state  4  occurs  after  all  degrees  of  freedom  have  relaxed 
to  their  equilibrium  values.  The  law  of  mass  action  for  the  second  reac¬ 
tion  of  Equation  (1)  is  Equation  (8) .  This  equation  together  with  Equa¬ 
tions  (7) ,  (15) ,  (16) ,  (17) ,  (11)  and  (13)  are  again  solved  by  iteration 
to  obtain  conditions  for  state  4. 


C  = 


T2  -  TX 
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For  the  fully  dispersed  wave,  state  1  and  2  are  equivalent;  that 
is,  there  is  no  normal  "shock  jump".  Since  the  initial  temperature  Tj 
was  low  in  our  experiments  and  since  fully  dispersed  waves  are  weak,  the 
second  reaction  was  not  excited,  and  states  3  and  4  therefore  coincided. 

In  order  to  describe  the  non-equilibrium  flow  field  in  detail,  it 
is  necessary  to  knew  the  time  dependent  behavior  of  the  nan-equilibrium 
variables  of  state.  These  relationships  are  given  by  the  rate  equations. 
Both  reactions  of  Equation  (1)  have  been  previously  investigated  and  the 
reaction  mechanisms  appear  to  be  known.  These  two  rate  equations  are  dis¬ 
cussed  in  Section  VI.  Utilizing  the  thermal  equation  of  state  (11) ,  with 
the  relationships  of  Equation  (6) ,  the  two  coupled  rate  equations  (61)  and 
(68)  can  be  rewritten  as. 


and 


-  4X  (a  -  8) 2 

dt  =  ^1  rF  [  (1  "  a)  "  kT  (1  +  aX  +  8X J  ] 

pi  o  o 


dt  ^20  kF  (a  ~  +  ltD2b 


2X  (a  -  B)2 

£ - 2 - + 

RT  1  +  otX  +  ex^ 
o  o 


da 

dt 


(23) 


(24) 


The  one  dimensional  steady  state  conservation  equations,  differ¬ 
entiated  with  respect  to  time  in  a  shock  fixed  coordinate  system  are. 


Mass: 

o 

li 

SI'S 

a 

(25) 

Momentum: 

dp  ,  M  dU 

5t  +  pUdt"  0 

(26) 

Energy: 

SIS' 

+ 

SN8 

II 

o 

(27) 

It  is  noted,  that  conditions  are  fixed  at  a  given  time  or  distance  down¬ 
stream  of  the  normal  shock  junp,  and  the  time  dependence  of  the  variables 
is  a  function  of  position  from  the  shock  front. 

Differentiating  the  thermal  equations  of  state  (11)  logarithmically 

yields 

d£  dW_d£._®=  (28) 

p  W  p  T 
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where 


aw  _  xo(da  +  361 
w  ”  l  +  ox  +  sx 

O  O 

and 

_ 3>  ...  ,da  dg.  1  dp  1  dT 

p  dt  "  1  +  aX  +  gX  \it  dtJ  p  dt  “  T  dt 

o  o 


Fran  the  caloric  equation  of  state,  (13) , 


(29) 


(30) 


dh  _  ,3h,  dT  ,3h.  da  ,3h.  dB 

dt  l3T'a,8  dt  ^3a^T,B  dt  '36  T, a  dt 


where  the  partial  derivatives,  evaluated  from  Equation  (13)  are. 


(31) 


(~)  =  C 

'3T'a,B  p 


(32) 


,3hv  _  Xo 

3a'T, B  W(1  +  aX  +  BX  )2 
o  o 


(13846  -  13174  8Xq) 


(33) 


and 


f3h,  _  Xo _ 

l3B;T,a  W(1  +  aX  +  BX  ) 2 
o  o 


(27020  +  13174  aX) 
o 


Finally, 


(34) 


^  -  0  =  0  (35) 


where  x  is  distance  measured  behind  the  shock  front. 

In  the  fixed  shock  coordinates  used,  the  time  t  is  the  flow  time. 
This  is  the  time  for  a  particle  to  travel  frcm  the  shock  front  to  soma 
position  x  downstream.  In  the  laboratory,  however,  time  is  measured 
from  the  passage  of  the  shock  front  across  the  observation  window  till  the 
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passage  of  a  particle  a  distance  x  behind  the  shock  front.  This  time, 
tlab  '  is  ec3ual  to; 


fclab 


x 

Ui 


(36) 


where  Uj  is  the  shock  velocity. 

These  two  times  are  related  through  Equations  (35)  and  (36)  and 
Figure  2  graphically  demonstrates  the  difference  between  these  two  times. 
It  is  seen  that  t  >  t^y,  .  This  causes  a  time  "compression"  on  oscillo¬ 
scope  traces  since  such  observations  of  e.g.  pressure,  record  changes  of 
state  at  a  fixed  location  in  the  shock  tube. 

Equations  (23)  through  (27) ,  (30) ,  (31)  and  (35)  comprise  eight 
coupled  first  order  non-linear  equations  in  the  variables  a  ,  6  ,  p  ,  U  , 
p  ,  h  ,  T  ,  and  x  .  These  equations  were  numerically  integrated  on  an 
IBM  7090-7094  DCS  computer  using  Milne's  second-difference  method  with  a 
Hamming  Predictor-Corrector.  (References  28,  29) . 

The  initial  conditions  for  these  eight  equations  cane  from  the 
measured  values  of  Xq  ,  pi  ,  Tx  and  Uj  from  which  ai  »  3i  ,  Pi  and 
hi  are  calculated.  The  integration  is  carried  out  until  a  and  8  are 
within  one  percent  of  their  equilibrium  values  at  state  4.  Since  the  flow 
variables  asymptotically  approach  their  final  values,  the  relaxation  time 
literally  goes  to  infinity.  It  is  physically  more  meaningful  therefore  to 
ascribe  a  value  to  x  when  the  reaction  has  reached  a  certain  stage  near 
the  final  equilibrium  state.  For  this  work,  x  was  arbitrarily  taken  as 
the  time  when  a  and  B  have  reached  99%  of  their  final  values,  there¬ 
fore. 


x 

ex 


=  0.99  34 


and 


(37) 


TB  -  fcB  =  0.99  84 


For  a  fully  dispersed  wave,  the  initial  conditions  are  the  values 
at  state  1.  A  difficulty  arises  since  this  is  an  equilibrium  state  and 
the  right  hand  side  of  Equation  (23)  is  equal  to  zero.  Hence,  a  will 
not  progress  to  its  value  at  state  4.  A  computational  device  was  used  to 
overcome  this.  The  initial  conditions  were  perturbed  slightly  and  the  in¬ 
tegration  then  proceeded.  The  perturbation  consisted  of  adding  10“ 6  to  the 
initial  value  of  a  .  The  effect  of  this  perturbation  was  shown  to  be 
negligible  by  doubling  the  perturbation  in  the  calculation  with  no  appre¬ 
ciable  effect  noted  in  the  result. 
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As  discussed  by  Vincent!  and  Kruger  (12)  a  continuous  solution  for 
a  partly  dispersed  wave  does  not  exist  if  the  initial  conditions  are  chosen 
as  those  of  state  1.  Therefore,  a  discontinuity  in  the  flow  variables  re¬ 
sulting  in  the  conditions  at  state  2  provides  the  initial  conditions  for 
this  case  and  the  equations  may  then  be  integrated  to  the  final  conditions 
at  t  g  and  t  g  • 

Ihe  concentration  of  N02*  is, 


[N02] 


p  “o*0  -  6) 

RT  1  +  oX  +  ex 

o  o 


(38) 


Fran  this  relationship,  the  ratio  of  the  concentration  of  N02  at 
any  point  in  the  flow  field  to  the  initial  concentration  of  N02  is  obtained. 
These  values  are  compared  to  the  experimental  values  in  Section  VII.  Also, 
the  calculated  values  for  several  experiments  of  a  ,  8  ,  T/T^  ,  p/pj  and 
U/Uj  vs.  t  are  given. 


*We  recall  the  definition  of  concentration  given  by 
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Figures  1.1  and  1.2  in  Appendix  I  shew  a  sketch  and  a  photograph  of 
the  shock  tube,  chemical  mixing  station,  and  instrumentation.  The  shock 
tube  is  28  feet  long  overall,  with  an  inside  diameter  of  3  inches.  The 
driver  section  is  4  feet  long  and  the  entire  tube  is  constructed  of  stain¬ 
less  steel  to  permit  the  use  of  the  highly  corrosive  nitric  oxides.  A 
jacket  surrounds  the  tube  so  that  the  walls  may  be  heated  or  cooled  by  cir¬ 
culation  of  suitable  liquids.  Thermocouples  monitor  the  temperature  along 
the  tube  length.  Shock  speeds  are  measured  by  timing  the  shock  passage 
between  two  optical  trigger  stations.  These  stations  use  an  "absorption" 
trigger  for  weak  shocks  and  a  conventional  schlieren  type  trigger  for 
stronger  shock  waves.  Between  these  two  stations  is  the  light  absorption 
station  which  measures  the  time  dependent  concentration  of  N02 .  In  the 
mixing  station,  the  reactant  gases  and  the  inert  gas  are  introduced  into  a 
cylinder  in  which  they  are  thoroughly  mixed  prior  to  introduction  into  the 
shock  tube.  A  detailed  description  of  the  equipment  is  given  in  Appendix  I. 

To  determine  the  non-equilibrium  flew  through  a  shock  wave,  it  is 
necessary  to  find  the  initial  conditions  of  the  flew  variables  entering  the 
wave  and  the  time  dependent  history  of  the  non-equilibrium  relaxation  pro¬ 
cesses. 


A  gas  mixture  is  prepared  as  described  in  Appendix  I,  giving  the 
initial  reactant  mole  fraction  Xq  .  This  mixture  enters  the  shock  tube 
where  the  pressure  is  recorded  to  ±  0.1  torr  and  the  temperature  to  ±  0.1  °K. 
The  shock  velocity  is  measured  by  timing  the  shock  passage  between  two  trig¬ 
ger  stations  as  described  in  Appendix  I.  Using  this  technique  the  shock 
velocity  is  knewn  to  better  than  0.3%. '  The  initial  conditions  Xq,  pi ,  T], 
and  Ui  are  thus  determined. 

The  time  dependent  concentration  of  one  species  of  the  reaction, 

[N02] ,  was  furthermore  measured  photometrically  to  determine  the  non¬ 
equilibrium  mode.  The  absorption  spectra  of  N204  and  N02  are  given  by  Hall 
and  Blacet  (30)  and  it  is  found  that  at  4350  A,  N02  absorbs  strongly  while 
N204  is  transparent.  The  inert  gases  nitrogen  or  argon  are,  of  course,  also 
transparent  at  this  wave  length  as  is  NO.  Figure  3  is  a  plot  of  I/Io  vs. 
[N02]  with  I/I0  as  the  ratio  of  transmitted  to  incident  light.  A  linear 
relationship  between  the  concentration  and  the  logarithm  of  I/IG  is 
observed.  Lambert-Beer's  law  is  therefore  experimentally  validated  and  we 
can  write 


_I_  =  10-e[NO2]d 
Xo 

where  d  is  the  optical  path  length. 
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Fig.  3  Intensity  Ratio  versus  N02  Concentration 
e  =  log10  (Io/I/[N021  d 
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Over  the  large  temperature  range  investigated  in  this  work 
(280  <  T  <  2100  °K) ,  the  molar  extinction  coefficient,  e  ,  was  found  to 
have  a  small  temperature  dependence.  However,  within  the  small  temperature 
changes  induced  by  the  weak  shock  waves  utilized  in  investigating  flew 
fields  with  only  one  non-equilibrium  mode,  i.e.  the  dissociation  of  N2C>4 , 
the  value  of  e  was  taken  to  be  constant. 

Using  the  slope  of  the  line  in  Figure  3  and  Equation  (39)  ,  a  value 
of  £  is  obtained,  and 


e  =  132.5  £/mole-an 


(40) 


in  the  temperature  range  280  <  T  <  310  °K.  These  measurements  were  made 
directly  in  the  shock  tube,  with  I  and  IQ  measured  as  explained  in 
Appendix  I,  and  with  [N02]  calculated  from  Equation  (38) .  Since  all  reac¬ 
tions  are  very  rapid,  the  system  remains  at  equilibrium  during  gas  handling 
operations  making  the  measurements  valid  as  described.  This  value  is  in 
agreement  with  that  of  Wegener  (4) .  A  series  of  neutral  density  filters 
was  employed  to  investigate  the  possibility  of  photo-dissociation.  It  was 
found  that  by  successively  reducing  the  incident  light  intensity,  IQ  ,  by 
a  factor  of  two,  the  output  signal,  I  ,  also  varied  by  the  same  factor  so 
that  I/I0  remained  constant.  Since  photo-dissociation  would  have  varied 
I/IQ  as  a  function  of  IG  ,  the  constant  ratio  of  I/IQ  shewed  photo¬ 
dissociation  to  be  negligible. 

At  the  higher  temperature  interval  (T  >  1200  °K)  in  which  the  sec¬ 
ond  part  of  reaction  (1)  takes  place  the  temperature  dependence  of  e  could 
no  longer  be  neglected.  Also,  e  could  not  be  measured  directly  as  it  could 
at  rocm  temperature .  This  is  due  to  the  fact  that  N02  is  unstable  at  these 
temperatures .  Therefore,  an  indirect  approach,  based  on  the  kinetic  record., 
of  the  dissociation  was  employed  to  find  the  molar  extinction  coefficient. 

In  the  following  relation; 


V3  _ 

V1  Il/Io 


(41) 


Vi  and  V3  are  the  measured  voltages  at  states  1  and  3  respectively. 

~L\/10  is  the  transmitted  to  incident  light  ratio  at  state  1,  and  I3/I0 
may  then  be  determined.  State  3  corresponds  to  the  maximum  concentration 
of  N02  as  described  in  Section  VII. 

Frcm  the  analysis  of  Section  III,  [N02]  3  can  be  calculated  at  tem¬ 
perature  T3  .  Now  losing  Equation  (39)  a  value  of  e  is  achieved  for  a 
specific  temperature. 

Figure  4  is  a  plot  of  z  vs.  T  frcm  a  series  of  experiments 
using  a  variety  of  conditions  of  state  3.  The  lack  of  dependence  of  e  on 
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Fig.  4  Temperature  Dependence  of  the  Extinction  Coefficient  at  4350A 
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concentration  of  N02  is  evidence  that  Lambert-Beer 1 s  law  is  still  valid 
at  the  higher  temperatures.  The  values  of  Schott  and  Davidson  (31) ,  and 
Huffman  and  Davidson  (8)  agree  well  with  this  work.  The  results  shewn  in 
Figure  4  may  be  represented  by  the  relationship 


logic  e  =  -  1.113  x  lO-4  T  +  2.161 


(42) 


for  a  range  of  temperatures  300  *  T  *  2100  °K. 

This  method  of  obtaining  e  has  the  effect  of  forcing  the  theore- . 
tical  and  experimental  values  of  [N02]  to  coincide  at  state  3.  For  the  low 
temperature  results  of  the  first  reaction  of  Equation  (1)  this  was  not  the 
case  since  e  was  determined  independently  of  the  analysis.  In  Section 
VII,  it  will  be  seen  that  the  theoretical  and  experimental  values  for  con¬ 
ditions  at  state  2  compare  well  for  the  first  reaction.  This  would  indicate 
that  the  prediction  of  conditions  at  state  3  for  the  second  reaction  essen¬ 
tial  for  the  measurement  of  e (T)  may  also  be  assumed  to  be  reliable. 

In  the  photometric  measurement  of  the  N02  concentration  a  narrow 
parallel  beam  of  light  is  employed.  This  beam  passes  at  right  angles  to 
the  tube  wall  through  the  gas  mixture  and  impinges  upon  the  photomultiplier 
tube.  Ihe  accuracy  of  these  measurements  depends  upon  two  factors,  the 
electronic  response  time  of  the  photomultiplier  system,  and  the  width  of  the 
parallel  light  beam.  The  electronic  response  time  of  the  photomultiplier 
system  used  was  measured  to  be  less  than  0.5  ys  (see  Appendix  I)  and  there¬ 
fore  it  is  much  more  rapid  than  the  phenomena  that  concern  us  here. 

In  Section  I  the  concept  of  the  Damkohler  parameter  was  introduced 
as  the  ratio  of  the  relaxation  time  to  a  characteristic  flew  time.  Simi¬ 
larly,  a  Damkohler  parameter  may  be  defined  in  terms  of  the  ratio  of  a 
relaxation  length  to  a  flow  length.  The  appropriate  flew  length  for  the 
measurements  made  in  this  system  is  the  slit  width  of  the  light  absorption 
measuring  system.  The  slit  width  was  here  chosen  to  be  large  enough  so  that 
D  <<  1  (equilibrium)  was  valid  for  translational,  rotational,  and  vibrational 
relaxation.  On  the  other  hand  the  slit  width  was  chosen  to  be  small  com¬ 
pared  to  the  chemical  relaxation  length. 

The  effect  of  the  slit  width  on  such  measurements  is  discussed  by 
Stupochenko  (1) .  With  a  finite  width,  6  ,  the  signal,  V  ,  measured  at  any 
given  instant  is  an  average  of  the  actual  signal,  I  ,  over  the  slit  width 
so  that 


V 


rt  +  6/2Ui 


I  dt 


-  6/2U! 


If  the  signal,  I 


is  approximated  by  an  exponential. 


(43) 
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I  =  IQ  exp  (  -  t/x) 


(44) 


where.  t  is  a  chemical  relaxation  time,  then. 


2tI  U,  .  . 

»»  o  1  —t/x  .  ,  <S 

V  =  — - - e  w  sirih  — 

o  2Ujx 


(45) 


Hois  equation  may  be  expanded  in  a  series  to  yield 


V  -  Ic  (1  +  2i  £lxl  + . )  (46) 


Hie  effect  of  the  slit  width  will  be  negligible  if  the  following  condition 
holds: 


<52 

24  11^x2 


«  1 


(47) 


In  our  experiments  the  slit  width  was  0.6  itm  and  if  an  extreme 
example  of  the  results  to  be  shown  later  is  considered,  values  of 
Ui  =  356  m  sec-1  and  x  =  3  ys  may  be  assumed.  For  this  situation,  the 
term  on  the  left  hand  side  of  Equation  (47)  is  equal  to  0.013.  This  value 
may  be  considered  negligible  compared  to  one.  For  the  above  shock  speed, 
the  time  required  for  the  shock  front  to  pass  across  the  slit  is  equal  to 
1.7  ys.  At  320  °K,  the  vibrational  relaxation  time  for  N204  is  0.3  ys 
(32,33) .  The  Danikohler  parameter  for  vibrational  relaxation  is  therefore 
equal  to  0.18.  This  value  of  D  is  sufficiently  small  to  justify  the  as¬ 
sumption  of  vibrational  equilibrium  as  stated  at  the  beginning  of  the  dis¬ 
cussion  of  the  effect  of  the  slit  width.  The  more  rapidly  adjusted  degrees 
of  freedom  of  rotation  and  translation  will  achieve  equilibrium  in  an  even 
shorter  period  of  time. 
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SECTION  V 

EXPERIMENTAL  RESULTS 


The  nan-equilibrium  flew  fields  were  experimentally  investigated 
using  the  following  technique.  A  gas  mixture  with  a  prescribed  mole  frac¬ 
tion  of  N2O4  at  a  point  where  there  is  no  dissociation,  Xq  ,  was  prepared. 
This  mixture  was  admitted  to  the  shock  tube  at  a  measured  pi  and  Tj  . 

Fran  these  values,  the  frozen  sound  speed  afi  of  the  mixture  was  computed 
by  Equation  (3) .  A  shock  wave  was  propagated  through  this  gas  mixture  and 
its  velocity  Ui  was  measured,  giving  the  shock  Mach  number  Mg  =  U^/afi  . 

As  shown  in  Figure  5,  oscilloscope  traces  of  voltage  versus  time  are  obtained 
in  which  the  voltage  is  proportional  to  the  transmitted  light  intensity. 

From  these  traces,  and  from  the  previously  measured  value  of  Iq  ,  the  time 
dependent  N02  concentration  may  be  computed  using  Equation  (39) . 

In  light  of  the  objectives  of  this  work,  the  experiments  were  di¬ 
vided  into  two  groups:  weak  shock  waves  in  which  only  one  non-equilibrium 
degree  of  freedom  was  excited;  and  shock  waves  with  two  non-equilibrium 
degrees  of  freedom. 

Figure  5  shows  oscilloscope  traces  for  fully  dispersed  and  partly 
dispersed  shock  waves  propagating  through  a  gas  mixture  of  N2C>4,  N02,  and 
A.  In  all  these  experiments,  the  conditions  at  state  1  were  the  same,  and 
the  shock  strengths  were  varied  by  changing  the  pressure  in  the  driver  sec¬ 
tion  of  the  shock  tube.  In  these  experiments  the  voltage  which  is  directly 
proportional  to  light  intensity,  is  inverted  so  that  the  positive  direction 
corresponds  to  N02  concentration  or  to  a  ,  the  progress  variable  or  degree 
of  dissociation.  The  first  three  experiments  (506,  505,  and  504)  give  re¬ 
sults  for  fully  dispersed  shock  waves  and  it  is  believed  that  they  represent 
the  first  demonstration  of  these  waves  for  a  chemically  relaxing  system. 

The  continuous  transition  frem  equilibrium  state  1  to  the  final  equilibrium 
state  is  clearly  seen.  The  next  six  experiments  involve  partly  dispersed 
shock  waves.  These  waves  continuously  increase  in  strength  as  seen  in  the 
traces  for  experiments  503,  500,  498,  497,  494  and  496.  Partly  dispersed 
waves  are  characterized  by  a  normal  shock  jump  across  which  only  the  trans¬ 
lational,  rotational,  and  vibrational  degrees  of  freedom  equilibrate.  This 
shock  junp  is  followed  by  the  chemical  relaxation  of  the  first  reaction  of 
Equation  (1)  to  its  final  equilibrium  value.  If  the  shock  waves  are  looked 
at  in  order  of  decreasing  shock  strength,  it  is  seen  that  as  predicted  be¬ 
fore,  the  normal  shock  junp  decreases  in  proportion  to  the  total  change  of 
state.  At  waves  propagating  at  the  frozen  sound  speed,  afi  ,  the  normal 
shock  jump  goes  to  zero.  Utilizing  the  methods  of  Section  III,  conditions 
can  be  measured  at  state  1,  the  initial  conditions  and  computed  at  state  2, 
conditions  after  the  normal  shock  jump,  and  at  state  4,  conditions  at  the 
final  equilibrium  state.  These  values  are  presented  in  Table  II.  1  of 
Appendix  II  for  the  nine  experiments  of  Figure  5.  In  Figures  12,  13  and 
14  in  Section  VII,  the  calculated  concentration  ratios  with  respect  to  con¬ 
ditions  at  state  1  are  compared  to  the  experimental  concentrations  obtained 
using  Equation  (39) .  The  agreement  is  seen  to  be  good. 
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Fig.  5  Fully  Dispersed  and  Partly  Dispersed  Shock  Waves  in  N2O4  -  NO2  -  A  at  Same  Initial  Conditions 
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The  second  series  of  experiments  is  concerned  with  shock  waves 
strong  enough  to  raise  the  temperature  sufficiently  to  cause  the  dissocia¬ 
tion  of  N02  in  the  second  reaction  of  Equation  (1) .  In  these  experiments 
it  was  found  expedient  to  vary  the  initial  pressure,  px  ,  in  the  shock 
tube.  This  had  the  effect  of  keeping  the  pressure  fairly  constant  after 
the  normal  shock  jump.  Figure  6  shews  oscilloscope  traces  for  these  strong 
shock  waves  for  two  series  of  experiments  using  gas  mixtures  of  different 
Xq  .  As  reported  in  Huffman  and  Davidson  (8) ,  an  orange  light  glow  was  at 
first  observed  to  emanate  from  the  shock  tube  for  these  strong  shock  waves. 
Since  the  carrier  gas  itself  and  the  reacting  gas  were  not  nearly  hot  enough 
to  radiate  at  the  shock  strengths  observed,  and  since  no  explanation  of  this 
phenomenon  was  available,  work  was  performed  to  determine  the  origin  of  the 
glow. 


At  first,  experiments  carried  out  with  pure  argon,  nitrogen  and 
helium  emitting  the  reactants  also  showed  this  glow.  This  demonstrates  that 
the  glow  clearly  was  not  caused  by  the  reacting  gas.  Experiments  with 
stranger  shocks,  up  to  Mg  =  8  showed  a  color  shift  of  the  radiation  to 
the  blue  indicating  a  continuum  black  body  radiation.  Calculations  on  the 
heat  transfer  through  the  boundary  layers  though,  shewed  that  the  walls  of 
the  tube  were  not  heated  above  100  °C  and  therefore  did  not  radiate.  By 
inserting  a  carefully  cleaned  glass  tube  into  the  shock  tube,  the  possibi¬ 
lity  of  radiating  scale  or  grease  on  the  walls  was  also  eliminated.  Experi¬ 
ments  were  then  performed  in  which  the  vacuum  pump  was  cold  trapped,  and 
research  grade  gases  of  high  purity  were  used.  In  these  experiments  the 
radiation  was  greatly  reduced  and  the  radiation  effect  was  therefore  shewn 
to  be  caused  by  the  organic  impurities,  such  as  pump  oil,  present  in  the 
tube. 


o  m  t 

In  the  experiments  shewn  in  Figure  6,  the  4350  A  absorption  filter 
and  a  series  of  collimating  slits  were  placed  between  the  photomultiplier 
tube  and  the  shock  tube  window.  In  this  configuration  with  the  research 
grade  gases,  no  radiation  was  detectable  and  the  effect  was  eliminated. 

In  Figure  6,  voltage  is  again  shewn  as  in  Figure  5  and  in  experi¬ 
ments  577  to  567  the  effects  of  increasing  the  shock  strengths  are  seen. 

A  normal  shock  jump  occurs  across  which  the  translational,  rotational,  and 
vibrational  degrees  of  freedom,  and  the  first  chemical  reaction  equilibrate. 
Following  the  normal  shock  jump  the  second  chemical  reaction  is  observed. 

As  the  shock  strength  increases,  as  expected  the  rate  of  dissociation  of 
NO 2  is  seen  to  increase.  Therefore  N02  disappears  and  consequently  the 
light  absorption  decreases.  In  experiment  577,  the  oscillations  at  the  end 
of  the  trace  can  be  shewn  to  be  due  to  the  arrival  of  the  contact  surface. 
This  was  not  observed  in  the  other  experiments  which  were  taken  at  higher 
sweep  speeds.  Experiments  587,  590,  596  and  599  shew  a  similar  trend  of 
dissociation  for  a  gas  mixture  at  a  lewer  value  of  Xq  . 

In  Table  II. 2  of  Appendix  II,  a  number  of  parameters  are  shown  which 
were  calculated  for  these  nine  experiments.  We  recall  that  state  1  cor¬ 
responds  to  the  initial  measured  conditions.  At  state  2,  the  translational, 
rotational,  and  vibrational  degrees  of  freedom  are  in  equilibrium,  while  the 
chemical  reactions  are  frozen.  At  state  3,  two  sets  of  conditions  are  given: 
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the  conditions  at  which  the  second  reaction  is  frozen  while  the  first 
reaction  goes  to  equilibrium,  and  the  conditions  at  which  the  second  reac¬ 
tion  proceeds  while  the  first  reaction  goes  to  equilibrium.  This  second 
set  of  values  corresponds  to  the  maximum  NO2  concentration  observed. 
Finally,  state  4  is  the  conditions  for  final  equilibrium.  Figure  18  shews 
the  experimental  versus  the  theoretical  values  of  the  NO2  concentration 
ratio  versus  x  ,  the  distance  downstream  of  the  normal  shock  front. 
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SECTION  VI 
CHEMICAL  KINETICS 


At  this  point,  all  experimental  parameters  have  been  measured  or 
computed  and  this  section  is  concerned  with  a  stud/  of  the  chemical  reac¬ 
tions  of  Equation  (1)  as  inferred  from  the  results.  For  the  two  reactions, 
the  relationship  ta  «  tg  is  valid  as  will  be  demonstrated  in  Section 
VII  and  therefore,  the  two  reactions  can  be  uncoupled.  In  the  analysis  that 
follows,  3  is  considered  frozen  in  the  first  reaction  which  takes  place 
at  temperatures  below  400  °K.  In  turn  for  the  second  reaction,  at  tempera¬ 
tures  above  600  °K,  a  is  in  equilibrium. 

A  rate  equation  for  a  chemical  reaction  may  be  of  the  general  form 


<H>t  =  "b  f  (p’  T'  “> 


(48) 


where  the  derivative  on  the  left  hand  side  refers  to  the  fact  that  the 
reaction  is  forced  to  take  place  at  a  constant  temperature.  In  such  an 
environment  kp  represents  a  constant.  However,  if  the  reaction  does  not 
take  place  at  a  constant  temperature,  which  is  the  case  for  this  work,  the 
temperature  dependence  of  the  rate  "constant"  kD  must  be  considered.  If 
a  rate  equation  is  of  the  form  of  Equation  (48) ,  then  the  coefficient,  kp  , 
may  be  found  from  the  oscilloscope  traces  which  give  the  time  rate  of  change 
of  the  concentration  of  NO  2  with  respect  to  a  fixed  observation  point  in  the 
shock  tube.  The  following  procedure  was  used  to  determine  kpi  as  defined 
by  Equation  (23) .  It  is  seen  that  this  equation  is  furthermore  of  the  form 
of  Equation  (48) .  later  in  this  section.  Equation  (23)  will  also  be  shown 
to  be  the  rate  equation  for  the  first  reaction  of  Equation  (1) .  If  we  de¬ 
fine  a  parameter 


5  =  1-  I/I. 


(49) 


then: 


da  d£d£ 
dt  dt/  da 


(50) 


Laboratory  time  is  related  to  the  flow  time  by  Equations  (35)  and  (36)  to 
give 


dt  d*'at]ab  ul 


dt|  ^  dx/dt  U 


(51) 
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Si  =  c  —+  {— ) 

da  p  da  3a 


(59) 


Equations  (54)  to  (59)  form  a  set  of  six  first  order  differential  equa¬ 
tions  which  may  be  numerically  solved  to  give  values  of  d£/da  ,  £  ,  p  , 

T  ,  and  U  as  functions  of  a  .  Now  using  Equations  (52)  and  (48) ,  values 
for  ko  are  finally  obtained  as  a  function  of  temperature. 

For  the  second  reaction  involving  the  dissociation  of  N02,  the  rate 
equation  is  more  carplex  and  contains  terms  expressing  two  carpeting  mech¬ 
anisms,  as  seen  in  Equation  (24) .  These  two  mechanisms  are  a  unimolecular 
and  a  bimolecular  process.  Since  Tg  >>  ra  ,  a  remains  at  its  final 
equilibrium  value  throughout  the  second  reaction,  and  the  rate  equation  has 
the  form 


'It1!  "  *D2u  fl  (P'  T'  S)  +  *D2b  f*  T<  6)  <60) 


The  procedure  used  in  obtaining  these  two  rate  constants  is  as  follows: 

The  derivatives  in  Equations  (52)  and  (54)  to  (59)  are  rewritten  with  re¬ 
spect  to  3  instead  of  a  .  These  equations  are  again  solved  numerically 
to  give  values  of  d£/dB  ,  B  ,  p  ,  T  ,  and  U  ,as  functions  of  B  .  Two 
limiting  cases  may  now  be  calculated  from  Equation  (60) .  In  the  first  case, 
kQ2b  is  set  equal  to  zero  so  that  the  rate  equation  reduces  to  the  form  of 
Equation  (48) .  Now  using  the  S  equivalent  of  Equation  (52)  and  the  first 
term  On  the  right  hand  side  of  Equation  (60) ,  values  for  kj)2u  are  obtained. 
Similarly,  ko2u  may  be  set  equal  to  zero  as  the  other  limiting  case  and, 
in  turn,  values  of  ko2b  maY  be  obtained. 

The  physical  significance  of  the  results  obtained  by  these  methods 
will  be  discussed  next.  Both  of  the  chemical  reactions  in  Equation  (1)  have 
been  extensively  investigated  in  certain  ranges  and  the  reaction  mechanisms 
appear  to  be  well  known.  The  pioneering  work  on  the  first  reaction  was  done 
by  Carrington  and  Davidson  (3).  It  is  interesting  to  note  that  the  study  of 
this  reaction  represented  the  first  application  of  a  shock  tube  to  research 
in  chemical  kinetics.  The  dissociation  rate  constant  was  found  in  a  manner 
similar  to  that  employed  here.  Wegener  (2)  utilized  the  inverse  technique, 
i.e.  the  cooling  of  the  gases  in  a  supersonic  nozzle  to  measure  the  re¬ 
combination  rate  constant  in  the  first  application  of  supersonic  nozzles  to 
studies  of  chemical  kinetics.  Other  work  on  this  reaction  was  performed  by 
Bauer  and  Gustavson  (33)  who  measured  relaxation  times  in  a  Kantrcwitz  type 
inpact  tube,  and  by  Sessler  (34)  and  Cher  (5)  who  measured  the  dissociation 
rate  constant  using  the  classical  method  of  ultrasonic  absorption. 

The  different  techniques  mentioned  have  given  results  in  good  agree¬ 
ment  and  the  rate  equation  deduced  from  the  experiments  is  given  by, 

3[N204) 

( - —)T  =  -  kjj!  (N204][M]  +  k^  [N02]2[M]  (61) 
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Here,  the  brackets  denote  concentration  and  M  stands  for  any  molecule 
participating  in  the  reaction.  Furthermore,  kD^  and  k^i  denote  the  dis¬ 
sociation  and  recaribinatian  rate  constants  respectively.  At  equilibrivm, 
3[N20iJ/3t  =  0  ,  and  the  right  hand  side  of  Equation  (61)  reduces  to  the 
law  of  mass  action  (7)  so  that 


(62) 


Here,  Jqji  has  units  of  liters  mole-1  sec"1,  k^i  has  units  of  liters2 
mole-2  sec-1  so  that  the  equilibrium  constant  has  units  of  mole 

liter-1.  Kpi  has  units  of  atmospheres  so  that  the  universal  gas  constant 
is  R  =  0.08205  liter  atm.  mole-1  °K-1.  Reaction  (61)  exhibits  a  second 
order  rate  of  dissociation  and  a  third  order  recombination  rate.*  The  dis¬ 
sociation  rate  constant  was  measured  using  oscilloscope  traces  similar  to 
those  of  Figure  5.  Utilizing  the  measured  values  of  5  and  d^/dtiab  ob¬ 
tained  from  these  traces,  koi  was  calculated  by  the  methods  presented 
above. 


Figures  7  and  8  are  plots  of  k^i  vs.  1/T  for  inert  gases  of 
nitrogen  and  argon  respectively.  The  data  is  seen  to  fit  the  Arrhenius 
equation:**  For  N2 


kDl  =  2.85  x  10 14  exp  (-  11,000/RT)  £, /mole-sec 


and  for  A  ; 


k^  =  2.23  x  10 14  exp  (-11,000/RT)  it/mole- sec 


(63) 


(64) 


Here,  as  is  customary  R  =  1.987  cal  mole  1  °K-1  . 

■These  results,  in  accord  with  those  of  Cher  show  that  nitrogen  was 
more  efficient  in  causing  a  reaction  than  argon.  The  results  of  this  work 
are  seen  to  be  in  good  agreement  with  those  of  Carrington  and  Davidson  (3) , 
and  Cher  (5) .  Errplcying  Equation  (62)  reconbination  rate  constants  may  be 
calculated  from  the  measured  Jqj  and  the  known  .  These  results  are 


*We  recall  that  the  order  of  a  reaction  may  be  inferred  by  adding  the  ex¬ 
ponents  in  the  concentration  terms. 

**In  equations  such  as  (63)  and  (64) ,  the  exponential  constant  refers  to  the 
activation  energy.  The  pre-exponential  term  is  a  function  of  the  collision 
rate,  and  the  probability  that  a  sufficiently  energetic  collision  will 
actually  result  in  a  reaction.  This  last  term  is  also  known  as  the  steric 
factor  (35) . 


29 


10-7  £/mole-sec 


AEDC-TR-69-229 


Fjg.  7  Dissociation  Rate  Constant  versus  1/T  for 
N204  +  N2  ^  2NOz  +  N2 
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Fig.  8  Dissociation  Rate  Constant  versus  1/T  for 
IM2O4  +  A  ^  2N02  +  A 
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plotted  vs.  tenperature  in  Figure  9  and  agreement  with  the  measurements  of 
Wegener  (2)  is  found.  We  note  that  kjQ  is  plotted  on  a  li  npar  scale  em¬ 
phasizing  the  cus ternary  scatter  of  such  findings. 

The  original  work  on  the  second  reaction  of  Equation  (1)  was  done 
by  Bodenstein  and  Rams  tetter  (6)  using  a  furnace  heated  reaction  vessel  in 
the  temperature  range  592  to  656  °K.  Similar  experiments  were  carried  out 
by  Rosser  and  Wise  (7)  in  the  temperature  range  630  to  1020  °K,  and  by 
Ashmore  and  Burnett  (9)  in  the  temperature  range  473  to  707  °K.  Huffman 
and  Davidson  (8)  also  investigated  this  reaction  in  a  shock  tube  in  the 
temperature  range  1400  to  2300  °K.  Our  work  covered  a  temperature  range  of 
1300  to  2100  °K.  The  results  of  all  these  efforts  indicated  that  the  reac¬ 
tion  mechanism  for  this  reaction  was  not  as  straightforward  as  for  the  first 
reaction.  Bodenstein  and  Rams  tetter  (6) ,  and  Rosser  and  Wise  (7)  suggested 
a  simple  bimolecular  reaction, 


3[N02] 

*  5t  *T  =  "  (65) 


where  the  value  kp2b  =  4.0  x  10 9  exp  (-  26,900/RT)  S. /mole-sec  was  obtained. 
Ashmore  and  Burnett  (9)  observed  a  faster  initial  rate  of  dissociation  than 
the  previous  two  investigators.  But  as  the  reaction  proceeded  and  NO  was 
produced  they  achieved  a  similar  result.  The  anomalously  fast  rate  was 
ascribed  to  the  additional  following  mechanisms: 


k2 

2N02  N03  +  NO  (66) 

*3 

and 


ki* 

NO  3  +  N02  — ►  N02  +  NO  +  02 


(67) 


Huffman  and  Davidson  (8)  proposed  the  result 


<■ 


3[N02] 

3t 


-  ‘  *b2b  [N°2j2  -  kD2u  I^HM] 


(68) 


in  which  the  second  term  is  important  at  high  temperatures  and  low  N02  con¬ 
centrations  .  Values  of  the  rate  constants  of  kn2b  =  2.5  x  10 10  exp 
(-  25,000/RT)  £/mole-sec  and  ko2u  =  3.06  x  101*  exp  (-  65,400/RT) 
a /mole- sec  were  proposed.  This  value  of  kp2b  is  eight  timeg  higher  than 
the  previous  values,  hewever ,  no  explanation  of  the  discrepancy  was  given. 
In  this  work,  a  rate  equation  of  the  form  of  Equation  (68)  was  assumed,  but 
somewhat  different  values  of  the  rate  constants  were  obtained  /than  those  of 
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Huffman  and  Davidson  (8) . 

In  Equation  (68)  both  terms  on  the  right  hand  side  are  of  second 
order.  Following  Huffman  and  Davidson  (8)  they  may  also  be  denoted  as  the 
brmolecular  and  the  unimolecular  dissociation  terms  respectively.  For  the 
high  temperatures  observed  for  this  reaction,  the  effect  of  reocmbination 
was  negligible  and  it  could  be  neglected. 

Figure  10  includes  a  plot  of  kD2b  vs.  1/T  determined  frcm  our 
experiments  by  assuming  a  rate  equation  of  the  farm  preposed  by  Rosser  and 
Wise  (7) ,  Equation  (65) .  It  is  seen  that  these  points  approach  the  values 
of  Rosser  and  Wise  shown  by  a  curve  at  the  low  temperature  end  of  the  scale 
(large  values  of  1/T) ,  but  they  are  much  too  large  for  higher  temperatures - 

The  experimental  constant,  in  expressions  for  the  rate  constants 
such  as  (63)  and  (64) ,  refers  to  the  activation  energy  for  a  specific  reac¬ 
tion.  Accordingly,  the  results  of  Rosser  and  Wise  (7)  corresponds  to  the 
activation  energy  of  the  bimolecular  reaction.  Since  the  points  of  Figure 
10  indicate  a  much  larger  value  of  the  activation  energy  at  high  tempera¬ 
tures  (i.e.,  an  increased  value  of  the  slope)  a  second  reaction  was  con¬ 
sidered  corresponding  to  a  higher  value  of  the  activation  energy.  This 
second  reaction  path  would  dominate  the  bimolecular  reaction  at  high  tem¬ 
peratures,  but  would  be  insignificant  in  the  temperature  range  considered 
by  Rosser  and  Wise  (7) . 

We  note  that  in  Figure  10  and  in  Figure  11,  the  solid  1  i  nPs  repre¬ 
sent  the  values  of  the  rate  constants  arrived  at  by  a  trial  and  error  pro¬ 
cedure  discussed  later.  It  is  also  noted  that  the  points  in  Figures  10 
and  11  serve  to  indicate  the  form  of  the  rate  equation  rather  than  to  give 
values  of  the  rate  constants. 

Following  Huffman  and  Davidson,  the  mechanism  for  a  unimolecular 
reaction  was  examined  suggesting, 


and 


N02  +  M 


NO  +  0  +  M 


k5 

0  +  N02  — >  02  +  NO  (fast) 


(69) 


(70) 


Reaction  (69)  gives  the  unimolecular  dissociation  mechanism  while  reaction 
(70)  is  sufficiently  fast  (t  <  8  ys)  to  be  considered  in  equilibrium  (8) . 
Accordingly,  the  rate  equation  for  reaction  (69)  is  proposed  to  be, 

3  [NO  ) 

<— at- JT  =  _  ’'dZu  Eno231m1  (71) 

where  M  is  any  particle  participating  in  the  reaction.  Using  this  equation. 
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Fig.  10  Bi molecular  Rate  Constant  versus  1/T  for 

1  dtNOgl 
D2b  [NO2I2  dt 
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the  results  of  Figure  10  were  recompiled  and  they  are  shown  in  Figure  11. 
Here,  the  slope  of  the  line  drawn  through  the  high  temperature  points  cor¬ 
responds  to  the  activation  energy  predicted  for  the  unimolecular  process. 
The  experimental  points  therefore  suggest  that  the  rate  equation  for  the 
complete  reaction  is  indeed  that  of  Equation  (68) . 

As  mentioned,  the  rate  constants  themselves  were  arrived  at  by 
using  a  trial  and  error  procedure.  Initial  values  of  the  rate  constants 
were  chosen  as  the  as  symp totes  of  the  two  sets  of  points  in  Figures  10  and 
11.  Using  the  theory  presented  in  Section  III,  the  theoretical  N02  concen¬ 
tration  was  calculated  as  a  f •unction  of  t]^  .  These  values  were  then 
compared  directly  with  the  experimental  traces  of  Figure  6.  The  rate  con¬ 
stants  were  then  adjusted  until  the  calculated  results  fit  the  experimental 
traces  for  both  high  and  low  temperature  experiments.  The  details  of  this 
technique  are  presented  in  Section  VII.  The  rate  constants  thus  derived 
were, 


=  3.0  x  109  exp  (-  26,900/RT)  5,/mole-sec  (72) 

and 

=  3*60  x  1013  exp  (-  65,400/RT)  5,/mole-sec  (73) 


These  are  plotted  as  the  solid  lines  in  Figures  10  and  11  respectively. 
These  values  are  in  agreement  with  those  of  Rosser  and  Wise  (7)  for  the  lew 
temperature  limit,  and  with  those  of  Huffman  and  Davidson  (8)  for  the  high 
temperature  limit. 
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SECTION  VII 

NON-EQUILIBRIUM  FLOW  FIELDS  WITH  CHEMICAL  REACTIONS 


When  the  values  of  the  rate  constants  are  known  together  with  the 
appropriate  rate  equations,  the  methods  of  Section  III  can  be  employed  to 
solve  the  conservation  equations  for  the  time  dependent  variables  of  state. 
The  resulting  set  of  eight  simultaneous  non-linear  first  order  differential 
equations  may  be  solved  numerically.  With  these  results.  Equation  (38)  pro¬ 
vides  the  N02  ooncentration  which  can  be  plotted  vs.  the  distance  downstream 
of  the  shock  front.  The  experimental  traces  and  Equation  (39)  similarly  can 
provide  a  plot  of  [N02]  vs.  x  where  x  =  Ui  tiah  . 

Figures  12,  13  and  14  give  the  N02  concentration  normalized  by  the 
N02  concentration  at  state  1  as  a  function  of  distanae  x  for  the  nine 
experiments  shewn  previously  in  Figure  5.  Argon  was  present  as  the  inert 
gas.  The  experiments  shown  are  for  shock  strengths  applicable  to  the  first 
reaction  of  Equation  (1)  only.  The  chief  source  of  error  in  obtaining  the 
rate  constants  from  the  oscilloscope  traces  is  incurred  in  the  measurements 
of  slope  of  the  photomultiplier  signal  with  a  certain  amount  of  electronic 
noise  present.  Figures  12,  13  and  14  serve  as  a  check  on  the  validity  of 
the  rate  constant  measurements,  and  they  help  to  eliminate  the  uncertainty 
of  the  slope  measurements.  Not  only  do  these  graphs  corroborate  the  values 
of  the  rate  constants,  but  they  justify  the  analytical  methods  of  Section 
III.  Figures  13  and  14  demonstrate  that  the  calculated  conditions  after  the 
normal  shock  jump  are  within  a  few  percent  of  the  experimental  values. 
Similarly,  the  conditions  for  final  equilibrium  are  in  good  agreement. 

Figure  12  shows  the  conditions  for  three  fully  dispersed  shock  waves 
of  different  strengths.  The  fit  between  theoretical  and  experimental  values 
is  excellent,  demonstrating  why  the  use  of  these  waves  is  ideal  for  a  study 
of  chemical  kinetics  as  suggested  in  Section  I.  Similarly,  Figures  13  and 
14  shew  the  experimental  vs.  the  theoretical  N02  concentration  ratio  for 
six  partly  dispersed  waves.  The  agreement  is  again  seen  to  be  good.  In 
experiment  496  though,  the  experimental  points  fall  slightly  be  lew  the 
theoretical  curve.  This  is  due  to  the  fact  that  the  limit  of  resolution  of 
the  equipment  is  being  reached  and  there  appears  to  be  an  electronic  lag  in 
the  signal. 

Figures  15,  16,  and  17  show  the  calculated  properties  of  a  ,  8  , 
p/pi  ,  T/Tj  ,  and  U/Ui  vs.  time  in  the  shock  fixed  coordinate  system.  For 
these  experiments  the  solutions  were  integrated  from  the  appropriate  equa¬ 
tions  from  the  initial  conditions  of  t  =  0  at  the  shock  front  to  t  =  Ta 
where  Ta  is  defined  in  Equation  (37) .  The  initial  properties  for  these 
graphs  are  given  in  Table  II.  1. 

Figure  15  shews  the  continuous  change  in  state  variables  for  a  fully 
dispersed  wave.  Figures  16  and  17  shew  the  change  of  variables  after  the 
normal  shock  front  at  t  =  0  .  The  behavior  of  the  temperature  ratio  is  of 
interest  in  that  it  first  rises  and  then  falls  to  its  final  value.  This  is 
due  to  the  energy  added  by  the  convective  motion  of  the  slewing  gas  flew  as 
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Fig.  12  Theoretical  and  Experimental  N02  Concentration  Ratios 
versus  Distance  for  Fully  Dispersed  Shock  Waves 
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[NC^  =  0.827  x  10“3  noles/liter  Px  =  200  Torr  Tj.  =  280.6  °K 

X  =  0.0793 
o 


Fig.  13  Theoretical  and  Experimental  NOz  Concentration  Ratios 
versus  Distance  for  Weak  Partly  Dispersed  Waves 
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Fig.  14  Theoretical  and  Experimental  N02  Concentration  Ratios 
versus  Distance  for  Partly  Dispersed  Shock  Waves 
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opposed  to  the  energy  subtracted  by  the  endothermic  chemical  reaction. 

For  sufficiently  strong  shocks ,  the  endothermic  reaction  dominates  and  the 
temperature  decreases  manotonically. 

A  flow  field  with  more  than  one  non-equilibrium  mode  may  be  experi¬ 
mentally  investigated  by  one  of  several  methods.  If  the  relaxation  times 
for  these  degrees  of  freedom  are  of  the  same  order  of  magnitude,  and  the 
reactions  occur  in  the  same  temperature  interval,  then  it  is  necessary  to 
measure  a  time  dependent  variable  of  state  for  each  non-equilibrium  mode. 

For  exanple,  in  a  chemical  relaxing  flew  with  two  simultaneous  reactions  it 
would  be  necessary  to  concurrently  measure  two  state  variables  such  as  the 
concentrations  of  two  species  of  the  reaction.  In  this  work,  we  have  the 
advantage  that  the  two  reactions  can  be  uncoupled  since  tg  >>  Ta  ,  and  it 
was  only  necessary  to  measure  one  progress  variable  i.e.  the  time  dependent 
N02  concentration .  Because  of  the  large  temperature  changes  introduced  by 
strong  shock  waves,  it  was  found  necessary  as  described  to  account  for  the 
temperature  dependence  of  the  extinction  coefficient  e  in  Equation  (39) . 
Due  to  this  it  was  necessary  to  measure  the  rate  constants  and  the  extinc¬ 
tion  coefficient  simultaneously.  An  iterative  procedure  was  employed  in 
which  values  of  the  rate  constants  were  assumed  and  adjusted  by  ccmparing 
the  calculated  and  experimental  values  of  the  time  dependent  N02  concentra¬ 
tions.  Finally,  the  calculated  conditions  at  state  3  and  Equation  (39)  were 
used  to  obtain  the  values  of  e(T)  shown  in  Figure  4. 

Figure  18  shows  the  theoretical  and  experimental  N02  concentration 
ratios  for  four  experiments  utilizing  the  values  of  Equations  (42) ,  (64) , 
(72),  and  (73).  At  state  3,  values  of  [N02]  g3/[N02]  x  were  calculated  as¬ 
suming  that  the  second  reaction  is  frozen  while  the  first  reaction  goes  to 
equilibrium.  Since  the  initial  rate  of  dissociation  of  N02  is  fairly  rapid 
the  maximum  value  of  [N02]  3/[N02]  i  obtained,  as  computed  kinetically,  is 
less  than  the  frozen  value  (Figure  18) .  Therefore,  state  3  was  defined  as 
the  actual  conditions  at  which  the  first  reaction  goes  practically  to  equi¬ 
librium.  This  state  corresponded  to  the  maximum  N02  concentration.  Since 
the  first  reaction  reached  equilibrium  within  0.1  ys,  the  assumption  that 
state  3  exists  immediately  following  the  normal  shock  jump  was  justified. 

In  calculating  conditions  at  state  3,  the  Ashmore  and  Burnett  (9)  reaction 
mechanisms  of  Equations  (66)  and  (67)  were  not  considered.  This  was  justi¬ 
fied  since  the  mechanism  of  Equation  (69)  dominates  the  reaction  at  high 
temperatures .  In  Table  II. 2  a  comparison  of  the  values  at  state  3  for  6 
either  frozen  or  relaxed  is  given.  Finally,  Figure  19  shows  several  cal¬ 
culated  properties  for  experiment  599. 
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Fig.  18  Theoretical  and  Experimental  N02  Concentration  Ratios  versus  Distance 
for  Strong  Partly  Dispersed  Shock  Waves.  Conditions  Given  after  Normal 
Shock  for  Beta  Frozen  and  Relaxed 
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SECTION  VIII 
SUMMARY 


Ihis  work  is  concerned  with  the  non-equilibrium  flow  of  a  chemi¬ 
cally  reacting  gas  mixture  through  a  shock  wave  in  a  shock  tube.  Ttoo  chemi¬ 
cal  reactions:  N204  +  M  £  2N02  +  M£2N0  +  02+M,  were  observed  utilizing 
a  light  absorption  technique  permitting  the  time  dependent  concentration  of 
the  species  N02  to  be  determined.  Argon  and  nitrogen  were  used  as  the  inert 
carrier  gas  M  .  The  second  of  the  two  reactions  studied  has  a  relaxation 
time  several  orders  of  magnitude  larger  than  the  first  reaction  and  moreover 
it  takes  place  at  higher  temperatures,  ftiis  allowed  the  two  reactions  to  be 
experimentally  uncoupled.  In  turn,  both  chemical  relaxation  Hires  were 
large  compared  to  the  translational,  rotational  and  vibrational  relaxation 
times.  Therefore,  all  molecular  degrees  of  freedom  other  than  those  per¬ 
taining  to  the  chemical  reactions  could  be  considered  in  equilibrium  through¬ 
out  the  flow  field. 

One  aim  of  this  work  was  the  investigation  of  the  passage  of  shock 
waves  in  a  shock  tube  through  a  gas  mixture  in  which  only  one  non-equilibrium 
degree  of  freedom  was  excited,  i.e. ,  the  first  of  the  above  chemical  reac¬ 
tions.  The  interplay  between  the  normal  shock  jump  and  the  chemical  relaxa¬ 
tion  zone  was  examined  as  a  function  of  shock  strengths.  This  produced  a 
complete  picture  of  the  flew  field  as  the  shock  strengths  were  varied  from 
weak  fully  dispersed  waves,  (propagating  be  lew  the  frozen  sound  speed) ,  to 
the  stronger  partly  dispersed  waves,  (propagating  above  the  frozen  sound 
speed) .  Fully  dispersed  waves  in  a  chemically  reacting  system  were  observed 
for  the  first  time  and  these  waves  were  found  to  be  ideally  suited  for  chemi¬ 
cal  kinetic  studies  due  to  the  inherently  long  relaxation  lengths  and  smal  1 
gradients  in  the  state  variables. 

For  stronger  shock  strengths,  the  second  reaction  was  observed  where 
the  two  reactions  were  uncoupled  due  to  the  large  difference  in  their  re¬ 
spective  relaxation  times.  Utilizing  this  technique,  a  composite  picture  of 
a  flew  field  with  two  non-equilibrium  degrees  of  freedom  was  formed.  Because 
these  two  reactions  could  be  examined  independently,  it  was  only  necessary 
to  measure  one  time  dependent  variable  of  state,  the  N02  concentration ,  in 
order  to  describe  the  flow  field. 

Finally,  a  study  of  the  reaction  mechanisms  and  rate  constants  for 
both  reactions  was  carried  out.  For  the  first  reaction,  following  Wegener 
(2)  the  rate  equation  was  taken  to  be  d[N204]/dt  =  -  k^  [N204]  [M]  +  k^ 

[N02 ] 2 [M]  .  For  argon  as  the  inert  gas,  the  rate  constant  was  found  to  be 
kjji  =  2.23  x  10 14  exp  (-  11000/RT)  Z  mole-1  sec-1  where  km  is  related  to 
koi  by  the  equilibrium  constant.  This  value  was  in  good  agreement  with  the 
values  of  other  investigators.  For  the  second  reaction,  following  Huffman 
and  Davidson  (8)  the  rate  equation  was  taken  to  be  d [N02  ] /dt  =  -  ko2b 
[N02  ] 2  -  k£)2u  [N02][M]  .  For  argon  again  taken  as  the  inert  gas  the  two 
dissociation  rate  constants  were  found  to  be  Jqj2b  =  3.0  x  10 9  exp  (-  26900/RT) 
Z  mole-1  sec-1  and  k^2u  =  3.60  x  1013  exp  (-  65400/KT)  z  mole"1  sec-1  . 

The  first  term  was  found  to  be  important  for  temperatures  belcw  1200  °K  while 
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the  second  term  became  dominant  for  higher  temperatures. 

An  analytical  description  of  the  flow  fields  was  performed  with  the 
assumptions  made  that  effects  of  viscosity,  heat  conduction,  and  diffusion, 
were  negligible  outside  the  classical  shock  region.  Also,  the  species  of 
the  gas  mixture  were  individually  considered  to  be  thermally  perfect.  The 
steady  state  conservation  equations  were  then  solved  numerically  together 
with  the  equations  of  state  and  the  two  appropriate  rate  equations.  This 
determined  the  time  dependent  histories  of  the  state  variables  in  the  flew. 
The  N02  concentrations  calculated  in  this  fashion  agreed  well  with  the  ex¬ 
perimental  values  determined  by  the  light  absorption  technique. 

By  observing  the  behavior  of  a  particular  well  kncwn  set  of  chemi¬ 
cal  reactions  in  a  flow  with  shock  waves  of  varying  strength  the  generalized 
behavior  of  a  non-equilibrium  flow  field  with  shock  waves  was  evidenced. 
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APPENDIX  I 

DESCRIPTION  OF  EQUIPMENT 


Figures  1.1  and  1.2  are  a  sketch  and  a  photograph  respectively  of 
the  shock  tube,  chemical  mixing  station,  and  instrumentation.  The  overall 
length  of  the  shock  tube  is  26  feet  8  inches  and  the  tube  is  constructed  of 
seamless  stainless  steel  304  tubing.  The  I.D.  of  the  tube  is  3.020  inches 
and  the  O.D.  is  3.25  inches,  except  for  the  test  section  whose  O.D.  is  3.50 
inches.  The  driven  tube  is  divided  into  three  sections  as  shown.  These 
sections  are  bolted  together  across  1  inch  inserts  which  hold  teflon  ”0" 
rings.  The  inside  of  the  tube  was  honed  to  insure  a  uniform  cross-section 
and  smooth  walls.  Around  each  length  of  the  driven  tube,  a  concentric  4 
inch  steel  pipe  has  been  welded,  forming  a  jacket  which  may  be  cooled  or 
heated  by  an  appropriate  fluid.  A  foam  rubber  insulation  covers  the  jacket 
and  the  manifold  connected  to  the  jacket.  Fourteen  copper-  cons  tan  tan  ther¬ 
mocouples  at  20  inch  intervals  are  inserted  in  the  jacket  to  monitor  the 
tube  temperature .  At  the  end  wall  a  thermocouple  is  inserted  directly  into 
the  tube  to  measure  the  gas  temperature. 

Flats  were  milled  on  the  test  section  for  three  sets  of  0.5  inch 
thick  quartz  windows.  The  flats  reduce  the  wall  thickness  of  the  tu t*3  at 
the  windows  to  1/32  inch  so  that  the  windows  lie  fairly  flush  with  the  in¬ 
side  diameter.  The  spacing  between  the  windows  was  measured  to  be  8.136  on, 
and  this  is  the  optical  length  used  in  the  absorption  measurements.  The 
two  outer  pairs  of  windows  are  used  for  the  passage  of  light  beams  acting 
as  trigger  stations  for  timing  the  shock  speed  and  initiating  the  scope  trig¬ 
ger.  They  are  1.00  meters  apart.  The  center  windows  are  used  for  the  light 
absorption  station. 

The  4  foot  driver  section  is  mounted  on  rollers  so  that  the  tube  may 
be  easily  opened  to  insert  the  diaphragms.  (For  certain  experiments,  a  six 
foot  driver  section  was  used) .  The  driver  section  is  attached  to  the  rest 
of  the  tube  by  fast  acting  clamps.  A  housing  is  welded  on  the  driver  sec¬ 
tion  containing  a  spring  loaded  firing  pin  used  to  rupture  the  diaphragms . 

The  entire  shock  tube  is  rigidly  mounted  on  a  welded  frame  of  4  inch  I  beam 
and  4  inch  pipe.  The  frame  is  bolted  to  the  floor  of  the  laboratory. 

The  driver  and  driven  sections  are  connected  to  the  chemical  mixing 
station  and  to  vacuum  pumps  by  0.5  inch  O.D.  stainless  steel  tubing.  The 
driver  section  is  independent  of  the  driven  section  so  that  it  may  be 
evacuated  and  filled  simultaneously  with  the  driven  section.  The  system  may 
be  connected  to  one  of  the  two  vacuum  punps,  a  Nash-Hitor  roughing  pump  or 
a  Welsh  Duo-seal  floor  putrp.  By  attaching  a  diffusion  pump  with  a  cold 
trap  the  system  could  be  pumped  down  to  10  microns.  The  outer  two  windows 
of  the  test  section  as  shown  in  Figure  1.1  are  used  for  shock  detection.  As 
described  by  Eberstein  (37)  an  optical  triggering  system  was  found  advisable 
for  the  weak  shocks  and  relatively  high  pressures  (above  30  torr)  of  interest. 
The  detection  system  may  be  used  in  one  of  tsro  configurations.  The  first 
configuration  is  a  standard  schlieren  deflection  system.  A  narrow  slightly 
converging  beam  of  light  is  focused  just  outside  a  0.125  inch  diameter  fiber 
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optic  light  guide.  The  light  guide  leads  to  a  RCA  931  A  photomultiplier 
tube.  The  density  change  across  the  shock  deflects  the  beam  into  the  light 
guide  and  the  resulting  signal  is  fed  into  a  pulse  amplifier  giving  a  re¬ 
sponse  time  to  the  shock  of  less  than  0.5  ys.  This  technique  was  not  suc¬ 
cessful  for  the  extremely  weak  fully  dispersed  shocks  and  at  pressures  be¬ 
low  80  ran  Hg  for  shocks  of  moderate  strength.  A  trigger  system  dependent 
on  the  absorptive  properties  of  N02  was  developed  for  these  situations .  In 
this  configuration,  a  4350  A  filter  is  placed  in  the  light  beam,  and  the 
beam  is  focused  directly  on  the  fiber  optic  light  guide.  Since  N02  strongly 
absorbs  light  at  this  wave  length,  the  increase  in  density  of  N02  through 
the  shock  causes  a  sharp  decrease  in  signal.  This  method  of  triggering  was 
successful  for  all  shock  strengths  and  pressures  used.  The  signal  from  the 
first  trigger  station  was  used  to  start:  the  type  5233L  Hewlett  Packard 
Electronic  Counter  and  to  trigger  the  delay  unit  in  the  oscilloscope.  The 
other  station  1  meter  downstream  stopped  the  counter  giving  the  shock  speed. 
The  absorption  station  as  shown  in  Figure  1.1  consists  of  a  narrow  parallel 
monochromatic  beam  of  light  passing  through  the  shock  tube  at  right  angles 
to  the  direction  of  fluid  flow  and  impinging  on  a  photomultiplier  tube.  The 
light  source  used  was  a  PEK  high  pressure  mercury  arc  lamp  type  107  with  a 
storage  battery  power  supply  to  insure  a  steady  DC  light  signal.  The  light 
was  collimated  by  a  series  of  slits  to  produce  a  parallel  beam.  The  light 
was  then  passed  through  a  Bausch  and  Lcmb  4350  A  absorption  filter  with  a 
half  width  of  50  A.  After  passage  through  the  tube,  the  light  entered  a 
RCA  IP  21  photomultiplier  tube  powered  by  a  Hamtjer  high  voltage  power  supply 
type  N401.  An  emitter-follower  was  used  to  increase  the  electronic  response 
time  and  the  resulting  signal  was  recorded  on  a  Tektronix  type  647A  oscil¬ 
loscope  with  a  built  in  delay  unit,  and  a  Polaroid  camera  attachment. 

A  high  speed  chopper  wheel  may  be  placed  in  the  light  path  and  it 
is  used  for  calibration  purposes  prior  to  a  run.  The  chopper  is  run  by  an 
air  driven  drill  turbine  and  it  is  capable  of  chopping  at  a  rate  of  20  Khz. 
Besides  being  used  to  measure  the  incident  DC  light  signal,  the  chopper  is 
used  in  conjunction  with  measurements  to  measure  either  the  electronic  re-  •' 
sponse  tirre  or  the  slit  width.  Since  the  measured  rise  time  is  a  ocmbination 
of  electronic  response  and  the  time  necessary  for  the  chopper  blade  to  pass 
through  the  light  beam,  by  measuring  the  rise  time  of  the  signal  for  a  known 
chopping  rate  the  electronic  rise  time  was  computed  for  a  known  slit  width. 
The  electronic  rise  time  was  of  the  order  of  0.5  us  which  was  sufficiently 
fast  compared  to  the  chemical  processes  studied.  Calibrated  neutral  density 
filters  may  also  be  placed  in  the  light  path  to  measure  the  linearity  of  the 
photomultiplier,  and  to  control  the  light  input. 

The  mixing  system,  as  shown  in  Figure  1.1  consists  of  a  tank  of 
volume  34.4  liters,  a  circulating  punp,  and  a  system  of  valves  to  admit 
gases  to  the  tank,  connect  vacuum  pumps  to  the  tank,  and  allow  the  gas  mix¬ 
ture  in  the  tanks  to  flow  into  the  shock  tube.  Due  to  the  poisonous  and 
corrosive  nature  of  the  gases  used,  the  components  of  the  mixing  station 
were  constructed  of  stainless  steel  and  teflon,  and  it  was  placed  in  a  walk- 
in  type  fume  hood.  The  temperature  in  the  tank  was  measured  with  a  copper- 
cons  tan  tan  thermocouple,  and  the  pressure  was  measured  roughly  with  a  bourdon 
type  gage  and  more  accurately  with  a  differential  stainless  steel  diaphragm 
strain  gage  pressure  transducer  (Dynisco  model  DPT  85-2) .  The  transducer 
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acts  as  a  null  gage  in  which  the  pressure  in  the  tanks  was  balanced  by  a 
manostat  and  the  pressure  was  read  on  a  Wallace  and  Tieman  type  FA  187 
mercury  manometer  or  on  a  Wallace  and  Tieman  pressure  gage  irodel  FA  145  to 
an  accuracy  of  0.1  torr.  This  technique  was  necessary  since  the  reacting 
gas  attacks  mercury. 

In  practice,  the  mixing  tank  was  evacuated  to  less  than  1  torr  and 
flushed  with  the  inert  gas  to  be  used  in  the  mixture.  The  procedure  was 
repeated  twice.  The  tank  was  then  filled  to  a  pressure  of  about  100  torr 
with  the  inert  gas.  The  N2O4  -  NO 2  was  then  admitted  to  achieve  the  de¬ 
sired  mole  fraction.  Inert  gas  was  then  added  to  a  pressure  of  2280  torr. 

At  each  stage  the  pressure  and  temperature  were  recorded  after  the  system 
had  been  allowed  to  reach  equilibrium.  The  gas  mixture  was  continually  cir¬ 
culated  through  the  tank  and  connecting  tubing  by  means  of  a  Randolf  purrp. 
Hus  ensured  that  the  mixture  was  homogeneous.  The  gases  used  were  supplied 
by  the  Matheson  Co.  Prepurified  nitrogen  with  a  minimum  purity  of  99.997% 
by  weight  and  high  purity  argon  of  purity  99.995%  were  used.  As  a  check, 
research  grade  argon  with  a  maximum  impurity  of  1  ppm  was  used.  The 
N2O4  “  NO 2  had  a  minimum  purity  of  95.5%  with  the  following  impurities: 
watfr.°*1%  max* '  nitrosyl  chloride  0.08%  max.,  nonvolitile  (ash)  0.01%  max. 
Additional  water  vapor  was  removed  by  passing  the  nitrogen  dioxide  gas 
through  a  silica  gel  drier.  To  check  the  purity  of  the  N2C>4  -  N02,  a  series 
of  experiments  were  performed  in  which  oxygen  was  allowed  to  circulate  for 
sight  hours  through  the  liquid  N2O4  in  order  to  oxydize  impurities.  This 
procedure  did  not  affect  the  properties  of  the  gas. 

The  procedure  followed  in  running  the  shock  tube  with  a  reacting 
gas  mixture  was  as  follows:  An  N2O4  -  N02  inert  gas  mixture  was  prepared 
as  previously  discussed.  A  diaphragm  of  cellophane  or  Mylar  was  placed  be¬ 
tween  the  driver  and  driven  sections  and  the  tube  was  clamped  together. 
Either  one  or  taro  sheets  of  diaphragm  were  used  depending  upon  the  pressure 
ratio.  The  diaphragm  was  near  its  bursting  pressure  when  it  ruptured.  The 
shock  tube  was  then  evacuated  to  less  than  0.5  torr  and  was  flushed  with  the 
inert  gas  and  pumped  down  again.  The  reactant  gas  mixture  was  next  admitted 
to  the  tube  and  its  pressure  and  temperature  were  measured  with  the  differ¬ 
ential  pressure  transducer  and  the  copper-cons tantan  thermocouple.  The 
driver  section  was  evacuated  and  the  driver  gas,  argon,  nitrogen,  or  helium 
was  brought  up  to  the  required  pressure.  For  very  weak  shocks,  a  diaphragm 
with  a  bursting  pressure  of  only  250  torr  was  used.  Therefore,  care  was 
taken  so  that  this  pressure  difference  was  not  exceeded  when  the  sections  of 
the  tube  were  being  evacuated  and  filled.  The  chopper  wheel  was  inserted  in 
the  light  path  and  the  incident  light  intensity  was  calibrated.  The  chopper 
wheel  was  removed  from  the  light  path  and  the  shock  was  initiated  by  re¬ 
leasing  the  firing  pin.  After  an  experiment,  the  reactant  gas  was  flushed 
out  of  the  tube  through  the  fume  hood.  This  was  done  by  filling  the  tube 
to  2280  torr  with  nitrogen  and  venting  it  through  the  hood.  This  procedure 
was  repeated  three  times.  The  tube  was  then  evacuated  twice  before  it  was 
opened  to  the  air.  After  every  few  runs  a  rag  soaked  in  acetone  was  pulled 
through  the  tube  to  clean  it,  then  a  chamois  cloth  was  pulled  through  to  re¬ 
move  any  remaining  dust  or  lint. 
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APPENDIX  II 

TABULATED  VALUES  OF  THE  CALCULATED  CONDITIONS 
ACROSS  THE  SHOCK  WAVE 


Tables  II. 1  and  II. 2  give  the  calculated  conditions  of  the  flew 
variables  after  the  normal  shock  jump  (state  2) ,  after  the  first  reaction 
of  Equation  (1)  has  reached  equilibrium  (state  3) ,  and  at  the  final  equili¬ 
brium  state  (state  4).  Table  II. 1  refers  to  the  experiments  of  Figure  5, 
and  Table  II. 2  refers  to  Figure  6. 

The  initial  conditions  at  state  1  refer  to  the  measured  values  of 
Pi  ,  Ti  ,  and  Ui  from  which  04  ,  3i  ,  and  [N02]  1  are  calculated. 

In  Table  II. 2,  two  values  are  given  for  conditions  after  the  first 
reaction  has  reached  equilibrium.  The  first  condition  was  calculated 
holding  3  frozen.  The  second  condition  is  the  value  at  state  3  with  3 
at  its  actual  value. 
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TABLE  11-1 


Final  Equilibrium 
State  4 


Initial  Conditions 
State  1 


After  Normal  Shock 
State  2 


On 


Exp  506  Mg  —  0.921  —  0.0793 


a 

e 

P  Torr 
T  °K 
U  nv/sec 

[N021  moles/liter 


0.4735 
0.68  x  10“4 
200.0 
280.6 
272.2 

0.827  x  10“3 


F.vp  505  Mg  =  0.961  XQ  —  0.0793 

0.4735 
0.68  x  10-4 
200.0 
280.6 
283.8 

0.827  x  HT3 
Exp  504  Mg  =  0.989  X^  =  0.0793 


a 

6 

P  Torr 
T  °K 
U  m/sec 

[N02]  moles/liter 


0.5231 
0.9  x  10“4 

230.6 

285.7 
241.3 

0.103  x  10"2 


0.5551 
0.11  x  10-3 
252.0 
289.0 
233.4 

0.118  x  10-2 


a 


0.4735 


0.5774 
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TABLE  11-1  Continued 


-j 


Initial  Conditions 
State  1 


After  Normal  Shock 
State  2 


Final  Equilibrium 
State  4 


6 

P  Torr 
T  °K 
U  nv/sec 

[NO2]  iroles/liter 


0.68  x  10-4 
200.0 
280.6 
292.1 

0.827  x  10“  3 


Exp  503  M  =  1.02 

s 


X  =  0.0793 
o 


a 

3 

P  Torr 
T  °K 
U  iVsec 

[NO2]  noles/liter 


0.4735 
0.68  x  lO"4 
200.0 
280.6 
299.9 

0.827  x  10-3 


0.4735 
0.68  x  10"4 

209.2 

285.3 

291.4 

0.851  x  IQ"3 


Exp  500  Mg  =  1.09 


X  =  0.0793 
o 


a 

3 

P  Torr 
T  °K 


0.4735 
0.68  x  10~4 
200.0 
280.6 


0.4735 
0.68  x  10~4 
247.3 
303.5 


0.13  x  10“3 
267.8 

291.3 

228.4 

0.129  x  10-2 


0.5980 
0.14  x  10“3 
282.8 

293.4 
223.9 

0.140  x  10-2 


0.6559 
0.20  x  10-3 
327.9 
299.6 
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TABLE  11-1  Continued 


Initial  Conditions  After  Normal  Shock  Final  Equilibrium 
State  1  State  2  State  4 


U  m/sec 

321.8 

281.3 

212.5 

[N02]  moles/liter 

0.827  x  10“  3 

0.946  x  10-3 

0.173  x  10-2 

Exp  498  Ms  =  1.10 

li 

x° 

0.0794 

a 

0.4790 

0.4790 

0.6730 

3 

0.68  x  10"4 

0.68  x  10"4 

0.22  x  10"3 

P  Torr 

200.0 

255.2 

337.3 

T  °K 

281.0 

307.5 

301.3 

U  m/sec 

326.5 

279.9 

210.6 

[N02]  moles/liter 

0.837  x  10“3 

0.976  x  10“3 

0.1821  x  10"2 

Exp  497  Ms  =  1.12 

ll 

x° 

0.0794 

a 

0.4790 

0.4790 

0.6860 

6 

0.68  x  10“4 

0.68  x  10“4 

0.24  x  10-3 

P  Torr 

200,0 

264.8 

348.3 

T  °K 

281.0 

311.7 

302.8 

U  ny'sec 

331.6 

278.0 

208.5 

[N02]  moles/liter 

0.837  x  10“  3 

0.998  x  10“3 

0.191  x  10"2 

TABLE  11-1  Concluded 


vo 


Initial  Conditions 
State  1 


After  Normal  Shock 
State  2 


Final  Eaui librium 
State  4 


Exp  494  M  =  1.19  X  =  0.0794 

~  s  o 


a 

6 

P  Torr 
T  °K 
U  m/sec 

[N02]  moles/liter 


0.4790 
0.68  x  10"4 
200.2 
281.0 

348.3 

0.837  x  10"  3 


0.4790 
0.68  x  10"4 

296.1 

325.4 

272.2 

0.107  x  10"2 


Exp  496  Mg  =  1.21  XQ  =  0.0794 


a 

6 

P  Torr 
T  °K 
U  m/sec 

[N02]  moles/liter 


0.4790 
0.68  x  10-4 
200.2 
281.0 
356.1 

0.837  x  10"  3 


0.4790 
0.68  x  10“4 

311.5 
331.8 
270.0 

0.110  x  10"  2 


0.7288 
0.31  x  10~3 

385.2 
307.7 
201.9 

0.220  x  10"2 


0.7487 
0.35  x  10“3 

403.3 
310.2 
199.0 

0.234  x  10'2 
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TABLE  II-: 


After  Alpha  Relaxation  Final  Equilibrium 
Beta  Frozen  State  3  State  4 


1.000 

0.32  x  10“  3 
945.0 
1305 
315.8 

0.174  x  10-2 


1.000 

0.162 

986.0 

1276 

302.4 

0.154  x  10"2 


1.000 

0.962 

1035 

994.6 

235.6 

0.881  x  10“4 


1.000 

0.35  x  lO-3 
1087 
1533 
334.9 

0.171  x  10“2 


1.000 

0.141 

1104 

1476 

321.0 

0.153  x  10-2 


1.000 

0.989 

1164 

1225 

268.3 

0.241  x  10“4 


1.000 


1.000 


1.000 
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TABLE  11-2  Continued 


'  Initial  Conditions 
State  1 


After  Normal  Shock 
State  2 


After  Alpha  Relaxation 
Beta  Frozen  State  3 


Final  Equilibrium 
State  4 


B 

0.35  x  10-3 

0.35  x  10-3 

0.35  x  10“  3 

0.126 

a.  993 

P  Torr 

49.00 

1040 

1038 

1076 

1102 

T  °K 

294.6 

1650 

1639 

1588 

1336 

U  m/sec 

1269 

341.1 

343.1 

332.4 

283.0 

[N02]  moles/liter 

0.356  x  10"  3 

0.132  x  10"2 

0.152  x  10"2 

0.139  x  10-2 

0.128 

Exp  566 

M  = 
s 

4.23  X  = 

o 

0.0799 

a 

0.8846 

0.8846 

1.000 

1.000 

1.000 

B 

0.37  x  10"  3 

0.37  x  10“  3 

0.37  x  10“  3 

0.116 

0.996 

P  Torr 

40.00 

898 

896.3 

1457 

943.6 

T  °K 

294.2 

1755 

1748 

1699 

1454 

U  m/sec 

1318 

350.1 

352.3 

341.8 

299.0 

[N02]  moles/liter 

0.288  x  10"  3 

0.108  x  10“  2 

0.122  x  10“2 

0.111  x  10— 2 

0.601 

Exp  567 

M  = 
s 

4.44  X  = 

o 

0.0799 

a 

0.8971 

0.8971 

1.000 

1.000 

1.000 

B 

0.40  x  10"3 

0.40  x  10"  3 

0.40  x  10“3 

0.103 

0.997 

\ 

P  Torr 

35.00 

869.5 

867.4 

875.0 

904.0 

T  °K 

294.3 

1906 

1901 

1855 

1617 
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TABLE  11-2  Continued 


Initial  Conditions 
State  1 


After  Normal  Shock 
State  2 


After  Alpha  Relaxation  Final  Equilibrium 
Beta  Frozen  State  3  State  4 


U  irv/sec 

1385 

361.0 

363.7 

354.6 

318.8 

[NOj]  moles/liter 

0.255  x  10-3 

0.978  x  10"3 

0.108  x  10“2 

0.995  x  10“  3 

0.315 

Exp  587 

Ms  = 

4.03  X  = 

0 

0.0431 

a 

0.9163 

0.9163 

1.000 

1.000 

1.000 

3 

0.44  x  10“  3 

0.44  x  lO-3 

0.44  x  10-3 

0.0841 

0.997 

P  Torr 

50.20 

1020 

1018 

1024 

1054 

T  °K 

294.6 

1687 

1683 

1663 

1516 

U  m/sec 

1274 

359.3 

360.0 

355.1 

326.2 

[NO2]  moles/liter 

0.208  x  10“  3 

0.736  x  10“3 

0.802  x  10“  3 

0.745  x  10“  3 

0.247 

Exp  590 

M  = 
s 

4.13  X  = 

0 

0.0431 

a 

0.9307 

0.9307 

1.000 

1.000 

1.000 

3 

0.48  x  10“  3 

0.48  x  10“  3 

0.48  x  10“  3 

0.0698 

0.998 

P  Torr 

40.30 

856.1 

856.2 

860.0 

883.5 

T  °K 

294.6 

1753 

1750.5 

1734 

1586 

U  m/sec 

1304 

365.0 

365.7 

361.7 

334.3 

[NO2]  moles/liter 

0.169  x  10“3 

0.604  x  10-3 

0.648  x  10-3 

0.610  x  10“  3 

0.148 
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TABLE  11-2  Concluded 


Initial  Conditions  After  Normal  Shock 
State  1  State  2 


After  Alpha  Relaxation 
Beta  Frozen  State  3 


Final  Equilibrium 
State  4 


Exp  596  =  4.31  Xo  =  0.0431 


a 

6 

P  Torr 
T  °K 
U  nv/sec 

w  [N02^  moles/liter 


0.9435 
0.54  x  10" 3 
34.90 

295.7 
1364 

0.148  x  10"3 


0.9435 
0.54  x  10“3 
810.0 
1893 

376.1 

0.536  x  10"3 


1.000 

0.54  x  10-3 

809.8 
1891 

376.9 

0.567  x  10”3 


1.000 

0.0570 

812.0 

1877 

373.8 

0.540  x  10“3 


1.000 

0.999 

830.8 

1731 

350.2 

0.859  x  lO"6 


Exp  599  Mg  =  4.44 


X  =  0.0431 
o 


a 

3 

P  Torr 
T  °k 
U  nv/sec 

[N02]  moles/liter 


0.9522 
0.58  x  10-3 
29.00 
295.7 
1406 

0.124  x  10“ 3 


0.9522 
0.58  x  10" 3 
716.0 
1993 

383.8 

0.454  x  lO"3 


1.000 

0.58  x  10" 3 
716.4 
1991 
384.6 

0.476  x  lO”3 


1.000 

0.0483 

717.1 

1979 

382.0 

0.456  x  10-3 


1.000 

0.999 

783.2 

1802 

351.1 

0.529  x  10-6 
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13  ABSTRACT 


The  present  work  is  concerned  with  a  theoretical  and  an  experiT 
mental  treatment  of  non-equilibrium  flow  of  a  chemically  reacting  gas 
mixture  through  a  shock  wave.  The  two  chemical  reactions  given  by 
N2O4  +  M  =  2NO2  +  M  =  2N0  +  O2  +  M,  where  M  represents  the  inert 
carrier  gases  argon  or  nitrogen,  were  studied.  The  experiments  were 
carried  out  in  a  temperature  controlled  shock  tube  and  a  light  absorp¬ 
tion  technique  permitted  the  time  dependent  concentration  of  the  species 
NO2  to  be  determined.  For  shock  strengths  where  the  temperature 
did  not  exceed  400°K  only  the  first  chemical  reaction  took  place. 
Stronger  shock  waves  excited  both  chemical  reactions  with  tempera¬ 
tures  up  to  2100°K,  and  for  this  reason,  flow  fields  with  two  non¬ 
equilibrium  modes  could  be  investigated.  Since  the  relaxation  times 
of  these  two  reactions  were  different  by  about  three  orders  of  magni¬ 
tude,  they  could  be  experimentally  uncoupled.  A  study  of  the  reaction 
mechanisms  and  rate  constants  for  both  reactions  was  carried  out. 

In  general ,  the  experimental  work  agreed  well  with  an  analytical 
description  of  the  non- equilibrium  flow  fields,  and  with  the  work  of 
previous  investigators  at  conditions  where  direct  comparisons  could 
be  made.  At  shock  strengths  exciting  only  the  first  chemical  reaction 
the  complete  picture  of  a  non-equilibrium  flow  field  with  only  one  non¬ 
equilibrium  mode  could  be  investigated.  In  this  situation,  the  shock 
strengths  were  varied  from  weak,  fully  dispersed  waves  to  strong, 
partly  dispersed  waves. 
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